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analysis  and  synthesis  of  speech  using  linear  predictive  coding  (LPC) . 

The  report  is  in  two  completely  separate  parts  which  were  originally 
the  reports  for  Masters  of  Science  projects  by  Paul  Hurst  (Part  I)  and 
Eugene  Evancoe  (Part  II) . 

Part  I  describes  an  approach  to  adaptive  autocorrelation  analysis  which 
could  be  implemented  using  analog  sampled  data  techniques  for  multi¬ 
plication  and  filtering. 

Part  II  is  an  investigation  of  the  use  of  high  order  switched- 
capacitor  ladder  filters  for  use  as  an  LPC  synthesis  filter. 
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1 .  Int 


Low  bit-rate  transmission  of  speech  has  been  an  important  speech 
research  goal.  Many  systems  have  been  realized  through  the  use  of  high 
speed  digital  computers.  Only  recently  has  the  aoplication  of  MOS-LSI 
made  possible  sophisticated  single  I.C.  speech  processors  [1],  At  this 
time,  the  most  visible  attempt  at  implementation  of  digital  transmission 
rate  reduction  are  the  "codec"  (coder-decoder)  chips.  Now  available  from 
a  number  of  companies,  the  codecs  use  a  logarithmic  transfer  curve  to 
achieve  a  digital  transmission  rate  of  6Hkbit3/sec .  This  is  a  large 
bandwidth  compared  to  the  bandwidth  of  speech  (4kHz).  Therefore,  it  is 
desirable  to  find  a  means  of  reducing  the  transmission  rate  even 
further. 


The  most  attractive  speech  bandwidth  compression  schemes  suffer 
from  the  need  of  fairly  complex  digital  processing  which  requires  expen¬ 
sive  and  complex  hardware.  In  this  project,  a  possible  solution  is 
explored  which  combines  analog  and  digital  circuitry  into  a  system  which 
can  operate  in  real  time  and  be  implemented  on  a  few  MOS-LSI  chips.  An 
approach  which  seems  promising  in  this  respect  is  an  analysis-svnthesis 
system  based  on  linear  prediction  (fig.  1).  Before  discussing  the  sys¬ 
tem,  a  short  review  of  linear  prediction  will  be  presented. 


*  A  similar  approach  is  applied  to  the  pitch  detection  problem  by 
Dalrymple  [21. 


Linear  Prediction  with  Infinite  Impulse  Response  Windows 
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Linear  prediction  is  applied  to  speech  to  extract  the  minimum 

information  to  describe  the  vocal  tract  frequency  response.*  Thi3  infor¬ 

mation  is  transmitted  digitally  at  a  very  low  rate  to  a  receiver  where 
the  speech  is  reconstructed  by  a  compatible  speech  synthesizer.  Such 
schemes  are  capable  of  transmitting  reasonaole  quality  speech  at  rates 
as  low  as  2. 4kb its/ second . 

Linear  predictive  coding  (L.P.C.)  of  speech  is  based  on  a  very  sim¬ 
ple  model  of  the  human  vocal  tract  (fig.  2  &  3).  Sounds  are  produced  by 
the  following  process: 

[1]  Air  is  forced  out  of  the  lungs. 

C2]  As  the  air  passes  through  the  glottis  (vocal  cords),  it  is  either 

allowed  to  flow  through  unrestricted  (unvoiced  sounds)  or  it  is 
transformed  into  short  bursts  of  air  by  the  periodic  opening  and 
closing  of  the  glottis  (voiced  sounds).  The  unrestricted  flow  of 
air  is  modelled  as  white  noise.  The  short  bursts  of  air  are 
modelled  as  impulses. 

[3]  The  air  then  passes  through  the  vocal  tract.  The  shape  of  the  vocal 
tract  and  the  position  of  the  lips  determine  the  resonant  frequen¬ 
cies  associated  with  the  sound.  (In  some  cases,  the  air  also  passes 
through  the  nasal  passage.) 

It  is  the  resonances  of  the  vocal  tract  that  we  attemDt  to  model 
using  linear  prediction.  The  shape  of  the  vocal  tract  varies  as 


*  An  excellent  introduction  to  linear  predictive  coding  is  given 
by  Makhoul  C  3 1  - 


Cross  Section  of  the  Vocal  Tract 
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different  sounds  are  produced,  but  It  varies  slowly  (every  10  to  20  mil¬ 
liseconds)  .  So  linear  orediction  can  be  applied  to  segments  of  speech  of 
this  length. 


L.P.C.  is  used  to  model  the  vocal  tract  resonances  by  fitting  an 
all-pole  transfer  function  to  the  vocal  tract  transfer  function  (fig. 
4).  The  model  has  the  form  (in  z-transform  notation): 


So<*> 

S^z) 


H(z)  = 


1-  l  a.z 
i=1 


-l 


p=number  of  poles  in  the  model . 


The  model  parameters  are  the  a^'s.  The  number  of  coles  in  the  model  can 
be  determined  by  the  following:  The  resonances  of  the  vocal  tract  are 
generally  separated  by  about  1kHz.  tach  resonance  requires  a  pair  of 
complex  conjugate  poles  in  the  model.  The  glottal  waveform  requires  1 
to  3  poles  for  modelling.  So  if  the  sampling  rate  is  f _ ,  the  model 
should  contain  at  least  p  poles,  where  p  is  given  by: 


P  > 


f 

s 

1kHz 


•4- 


1  . 


More  poles  improve  the  model  accuracy.  But  a  minimum  number  of  poles  i3 
desired  to  minimize  the  transmission  rate. 


One  approach  to  the  calculation  of  the  L.P.C.  coefficients  (the 
a^'s  in  the  vocal  tract  transfer  function)  is  the  autocorrelation 
approach  [3]>  Tnis  approach  requires  the  computation  of  p+1  autocorrela¬ 
tion  values  of  the  speech  waveform  computed  over  a  period  during  which 
the  speech  is  not  changing  (i.e.  only  one  sound  is  being  made).  The 
autocorrelation  values  can  be  directly  transformed  into  the  L.P.C.  coef¬ 
ficients  (see  Appendix  A).  To  accomplish  this,  the  speech  is  bandlimited 
to  half  the  sampling  frequency  and  then  sampled.  A  string  of 


"H\r  Vo  Ci>_l 
Cc<\^-S  ' 


F.wvx*  Ha) 


approx imatel y  240  consecutive  speech  samples  arc  then  used  in  the  auto¬ 
correlation  calculation: 

oo 

R(k)  =  l  s  (i)s  ( i*k)  k=0, 1 .... n 
i=-oo  w 

where  refers  to  th-  "windowed"  speecn  samples  (see  below).  The  model 
parameters  (the  a^'s)  can  be  calculated  directly  from  the  R(k)'s  (see 
appendix  A)  . 

The  usual  approach  is  to  multiply  the  .samples  of  the  incoming 
speech  waveform  by  a  finite  impulse  response  (F.I.R.)  time  window.  Typi¬ 
cal  windows  that  are  applied  to  speech  arc  Hanninoc  and  Hamming  windows 
[4],  A  muen  simpler  method  for  calculating  the  autocorrelation  values 
was  pointed  out  by  Barnwell  [5].  Rather  than  use  a  finite  imbulse 
response  time  window  on  the  sDeech,  Barnwell  considered  usinz  an  infin¬ 
ite  impulse  response  (I.I.R.)  time  window.  As  Ion?:  as  the  I.I.R.  window 
is  very  small  outside  a  30  msec  wide  rezion:  the  results  should  be  very 
similar  to  those  obtained  from  the  aDplication  of  a  F.I.R.  window.  The 
advantage  of  the  I.I.R.  approach  .is  that  the  I.I.R.  window  can  be  real¬ 
ized  by  a  simple  recursive  filter.  A  brief  summary  of  Barnwell's  work 

\ 

\ 

follows: 

Tne  speech  is  sampled.  Let's  call  the  sequence  of  sampled  soeech 
s(n).  This  input  sequence  is  divided  into  "frames"  consistinz  of 

ft 

240  consecutive  samples  (for  an  3kHz  sampling  rate).  Each  frame 
is  multiplied  by  a  time  window.  Let  the  index  "j"  refer  to  the 
largest  index  of  a  sample  in  any  frame  (i.e.  frame  j  ends  with  sam- 

*  Only  for  a  F.I.R.  window.  An  I.I.R.  window  multiplies  an  infin¬ 
ite  number  of  samples. 
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pie  a(j)  and  begins  with  33-1101®  s(i1-24'))).  For  convenience,  the 

ti-ne  window  w(l)  is  defined  such  that: 

w(i)*0  for  i<0 
w(i)=0  otherwise 


For  a  finite 

impulse  response  v  .dow, 

w(  i  Wo 

only 

for  -24CKi<0 

(fi«.  5). 

After  aDplyin^  the  time 

window 

to 

frame  j ,  a  new 

sequence  is  created: 


x(  i  ,.j)  =  s(  i)w(  j-i) .  (1) 

The  autocorrelation  values  for  frame  j  can  then  be  calculated  as 

00 

R(k,j)=  2  x(  i,j)x( i+k, j)  kzO,  1 ,2, . . .0.  {2} 

i=-oo 

where  H(k,j)  is  the  kth  autocorrelation  value  calculated  from  the 
samples  in  frame  .J. 

Usually  a  finite  length  Hannins  window  of  30  msec  width  is  used.  If 
a  time  window  which  is  infinite  in  length  but  very  small  outside  a 
30  msec  period  is  used  instead,  the  L.P.C.  system  is  greatly  sim¬ 
plified  . 


The  time  window  which  will  be  used  is  the  impulse  response  of  a 
second  order  filter  bavins?  two  coincident,  real  poles: 


W(z)  = 


n-cu-V  ‘ 


Its  time-reversed  impulse  response  is  (fi?.  5): 


w(n)=( 1-n)d-n  n^0 

w(n)=0  n>0 


Usin?  equations  (1)  and  {2},  we  qet 


* 
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jMV-*. 


oo 

rt(k,j)s  l  3( l)w( j-i) s( i+k)w( J-l-k) . 
i=-oo 

Now  define 


s'( i ,k)=s( i)s( i+k) 
w' ( i ,k}=w( i)w( i-k) . 


IV 


Then  we  can  write  equation  {3}  as 

00 

R(k,j)=  l  s'(i,k)w'( j-i,k) . 
i=-oo 

From  the  above  equation  it  can  be  seen  that  R(k,j)  is  the  convolu¬ 
tion  of  s'(i,k)  with  w'(i,k).  So  if  we  can  produce  the  sequence 
s'(i,k)  and  pass  it  through  the  linear,  time  invariant  filter  with 
impulse  response  w'(i,k),  we  will  have  calculated  the  autocorrela¬ 
tion  function  for  lag  k  (using  the  data  from  frame  ,j). 

Producing  the  sequence  s'(i,k)  requires  only  multiDl ication  and 
delay.  We  need  to  find  the  filter  corresponding  to  w'(i,k).  Recall 
that 


w'  ( i  ,k)=w( i) w( i+k) .  { 4 ) 

The  corresponding  z-transforms  are  defined  as: 

w(  i)  ->  W(  z) 
w(  i-k)  ->  z-l<W(z) 
w'  (  i  ,k)  W'  k(  z) 

Since  multiplication  in  the  time  domain  corresponds  to  convolution 
in  the  frequency,  domain,  we  can  get  W'  ( z)  from  equation  (U): 

W*  (z)  =  W(z)*z"kW(z)  =  J*  W(v)(— )"kW(— )v-1dv. 
k  2n  j  V  V  V 
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This  integral  cm  be  evaluated  to  give 


,k*2  -1 


V2>  - 


(I-*2.-1)3 


Note  that  each  value  of  k  corresponds  to  a  different  filter.  All 

2 

the  filters  have  3  Doles  at  d  and  one  real  zero.  Barnwell  has 
found  experimentally  that  d=0.<?8  is  the  best  choice  for  a  d  Dole 
L.P.C.  model  with  an  3kHz  sampling  rate. 


Notice  that  the  I.I.R.  system  does  not  require  any  lengthy  memory 
storage.  In  fact,  the  only  elements  of  the  system  are  a  sampler,  a  10- 
element  delay  line,  multipliers,  and  filters.  Rather  than  implement  the 
system  on  a  digital  computer,  it  is  possible  to  realize  the  entire 
L.P.C.  system  using  MOS  switched-capacitor  circuit  techniques  to  handle 
the  autocorrelation  computation  with  a  microprocessor  handling  Durbin's 
Recursion  (Appendix  A).  This  might  lead  to  a  3  or  U  chip  L.P.C.  speech 
analysis  system  (fig.  6). 

Before  designing  the  integrated  circuits,  some  important  questions 
ha/e  to  be  answered.  For  example:  How  good  must  the  multipliers  be?  How 
important  are  the  characteristics  of  the  10  integrating  filters?  What  is 
the  best  approach  to  the  system  which  will  perform  all  the  desired  com¬ 
putations  and  take  as  little  silicon  area  as  possible?  To  answer  these 
and  other  questions,  the  system  was  constructed  on  a  breadboard  U3ing 
only  standard  MOS  building  block3  (op  amps,  switches,  and  capacitors). 


r  n 

p\D  AC 

v  y 


r  \ 

pm  Ac. 

y  y 


A  breadooard  version  of  the  I.I.R.  system  was  built  using  the  10 
filters  suggested  by  Barnwell.  Following  his  results,  the  system 
operates  at  an  3kHz  sampling  rate  and  uses  a  9  pole  model. 

As  mentioned  earlier,  the  elements  of  the  system  are  1)  a  sample 
and  hold,  2)  a  delay  line,  3)  multipliers,  and  4)  filters.  Since  a 
breadboard  was  being  constructed,  circuits  employing  a  minimum  number  of 
precision  ratioed  capacitors  were  favored.  This  was  an  important  con¬ 
sideration  in  selecting  a  multiplier  scheme. 

All  the  individual  parts  of  the  system  will  be  presented 
separately.  The  following  discussion  assumes  that  all  circuit  elements 
are  ideal  and  that  all  capacitor  ratios  are  exactly  as  desired.  Non¬ 
idealities  will  be  considered  in  a  later  section. 

All  M03  transistors  drawn  in  the  figures  are  n-channel  devices. 
Therefore  they  act  like  a  closed  switch  when  their  gate  voltage  is  high. 
Also,  to  prevent  adjacent  switches  from  being  on  simultaneously,  all 
clock  signals  are  non-overlapping  (i.e.  and  are  non-overlapping  com¬ 
plements)  . 

Jl.J..  Sample  and  Hold 

A  sample  and  hold  consists  of  a  switch  for  sampling,  a  capacitor 
for  storing  the  sampled  voltage,  and  an  op  amp  for  buffering  (fig.  7). 
The  sample  and  hold  circuit  acquires  a  new  speech  sample  every  125 


microseconds. 
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3.-2*  Delay  Line 


This  particular  delay  line  uses  a  g-p'nase  clock  (fig. 7).  On  SAM¬ 
PLE:  1,  the  input  voltage  is  stored  on  C.,.  On  RESET: 1,  C2  discharges  com¬ 
pletely.  On  TRANSFERsI,  the  charge  on  C1  i3  forced  onto  C2-  This  gives 
an  output  voltage 


V  (n+l)  =  -  -r1  *  V.  { n) . 
out  in 


This  voltage  remains  at  the  output  during  the  next  SAMPLE  pulse,  allow¬ 
ing  similar  stages  to  be  cascaded  to  form  a  delay  line.  (Ideally,  C.,=C2 
which  gives  VQUt(n+1):-Vin(n) .) 


The  Multiplier 


A  high-accuracy  multiplier  was  desired  for  use  in  this  system.  Ana¬ 
log  multipliers  typically  have  accuracies  on  the  order  of  V4.  Rather 
than  build  a  breadboard  with  such  a  limited  multiplier  accuracy,  it  was 
decided  to  construct  the  multiplier  from  an  analog-to-digital  converter 
(A/D)  and  a  multiplying  digital-to-analog  converter  (MDAC).  This 
approach  can  provide  sufficient  accuracy  and  allows  flexibility  for 
experimentation  (i.e  a  simple  wiring  change  allows  us  to  see  how  the 

ft 

system  will  perform  with  an  3-bit  multiplier). 


*  Since  the  output  of  each  multiplier  is  oassed  through  a  low  pass 
filter  (bandwidth  of  about  30Hz)  ,  it  is  not  clear  that  a  high- 
accuracy  multiplier  is  required.  But  a  large  dynamic  range  is 
necessary. 


The  _A/D  Converter 

The  A/D  converter  used  was  first  presented  by  R.  McCharles  [6].  It 
requires  only  5  precision  ratioed  capacitors,  2  op  amps,  switches,  and 
some  associated  control  logic  (fig.  3).  The  digital  word  is  output  seri¬ 
ally.  The  digital  output  is  always  monotonic.  The  number  of  bits  per 
conversion  is  set  by  the  control  logic,  but  there  are  many  circuit  limi¬ 
tations  which  set  a  practical  upper  limit  on  the  number  of  cits. 

An  11-bit,  bipolar  converter  was  constructed  (10  bits  ♦  sign  bit). 
It  accepts  inputs  between  +10V  and  -10V.  The  digital  word  is  an  offset 
binary  representation  of  the  analog  voltage. 

The  McCharles'  A/D  Converter  implements  a  non-restoring  divide 
algorithm.  If  we  consider  each  op  amp  and  its  integrating  capacitor  a3 
an  analog  register,  the  operation  of  the  converter  can  be  written  com¬ 
pactly: 

1.  Clear  (Op  Amp  1)  and  (Op  Amp  2)  {discharge  all 

capacitors) 

2.  (Op  Amp  1)  4-  Vin1;  S=sign(0p  Amp  1)  {load  Op  Amp  1, 

set  sign  bit) 

3.  (Op  Amp  2)  4-  (Op  Amp  1)  {shift) 

4.  (Op  Amp  1)  4-  2*( Op  Amp  2)  -  S*Vref ;3=sign(  Op  Amp  1) 

{set  MSB) 

5.  for  i=1  to  9 

6.  (Op  Amp  2)  4-  (On  Amp  1)  {shift) 

7.  (Op  Amp  1)  4-  2*( Op  Amp  2)  -  S*Vr<lf  ;S=sign(0p  Amp  1) 

{set  next  bit) 

8.  next  i 

9.  end  (done  after  LSB) 

The  function  "sign(0p  amp  1)"  is  defined  to  be  the  logical  "0"  when 
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the  op  amp  output  voltage  Is  negative  arid  the  logical  ”1"  when  the  op 
amp  output  voltage  is  non-negat i ve .  The  variable  "S"  is  th®  serial  out¬ 
put  of  tne  A/D  converter.  One  conversion  requires  12  clock  cycles.  To 
complete  one  conversion  every  125  microseconds  (3kHz  sampling  rate),  the 
A/D  converter  uses  a  96kHz  clock. 


2  The  Multiplying  JO/A.  Converter 


A  simple  0/A  converter  was  designed  oy  H.  Fellman.  It  requires  only 
2  precision  ratioed  capacitors,  making  it  an  ideal  candidate  for  the 
breadboard.  It  is  capable  of  doin.^  four  quadrant  multiplication  and  is 
compatible  with  the  McCharles*  A/D  (fie?.  9). 


The  two  equal  capacitors  (C^  =  C2)  are  used  to  divide  the  voltage 
V|n2  by  successive  factors  of  two.  These  fractions  of  Vin2  are  then 
selectively  integrated  by  C,  depending  on  the  bit  of  the  incoming  digi¬ 
tal  word  (see  fig.  9  for  an  example).  After  12  clock  cycles,  is 
cnar<ted  to  a  voltage  equal  to 


V  .  ( n) 
out 


(n)'Vln2(n) 
’  10.0V 


2.4..  Filters 


The  system  requires  the  use  of  filters  at  both  the  inout  and  the 
output.  The  incoming  speech  must  be  band  limited  to  4kHz  and  pre¬ 
emphasized  before  it  is  processed.  The  output  of  the  multipliers  must  be 
integrated  by  the  10  low  pass  filters  to  give  the  autocorrelation 


values . 


i  LHI 


To  bandlimit  the  incoming  speech  to  4kHz,  a  fifth  order  Chebyshev 
filter  was  used.  It  has  Idb  of  ripple  in  the  passband  (0  to  3kHz).  At 
4kHz,  the  response  is  down  24db  (fi*.  10).  A  complete  discussion  of  the 
filter  can  be  found  in  (71. 

The  soeecn  was  pre-emphasized  to  whiten  the  spectrum.  This  results 
in  a  better  L.P.C.  model  [31.  The  pre-emoha3is  filter  has  one  zero  at 
500Hz  and  a  pole  at  4kHz  (fi<*.  11). 

2.- 4,.<£.  Outout  filters 

The  10  filters  at  the  output  of  the  multipliers  have  3  coincident, 
real  polos  and  one  real  zero.  Those  filters  can  be  realized  simply  (fiz. 
12).  The  first  two  sections  of  the  filter  each  realize  a  pole  at 
z=d^=0 . MoQU .  The  third  section  realizes  a  pole  and  a  zero.  The  location 
of  the  zero  depends  on  the  autocorrelation  laz  k.  Table  1  lists  the 
pole-zero  locations  for  all  10  filters.  As  with  3ll  switched-capacitor 
filters,  the  locations  of  tne  poles  and  zeros  are  determined  by  capacitor 
ratios. 

The  switched-capacitor  fiLters  all  nave  a  zain  of  1.00  at  D.C. 
Table  1  shows  that  all  th.e  filters  should  have  a  slizhtly  different  zain 
at  D.C.  This  D.C.  jain  factor  can  be  included  in  the  first  step  of  the 
Durbin's  Recursion  Algorithm  (see  appendix  A). 

The  low-pass  action  of  the  output  filters  reduces  the  contribution 
of  quantization  noise  and  circuit  noise  at  the  outputs. 
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4..  Circuit  Limitations 

On  the  breadboard ,  the  parts  used  were: 

silver-dipped  mica  capacitors 
glass  capacitors 
National  LF356  op  amps 
Siliconix  DG201  CMOS  switches 
Siliconix  DG190  JFET  switches 
Siliconix  DG39?  CMOS  switches 
Siliconix  D130  switch  drivers 
Signetics  SD5003  DMOS  switches 
National  LF393  sample  and  holds 
74LS00  series  logic  gates 

The  LF356  op  amps  have  a  JFET  input  stage,  so  the  input  bias 
currents  are  very  small  (<  8  nanoamps).  This  is  necessary  for  switched 
capacitor  circuits  since  the  input  bias  currant  must  flow  through  the 
capacitors.  These  op  amps  are  unconditionally  stable  when  driving  capa¬ 
citive  loads  of  10  mnofarads,  which  allows  large  capacitors  to  be  used. 
Other  excellent  properties  of  the  356  are  high  slew  rate,  low  noise,  and 
fast  settling. 

The  commercial  switches  all  have  very  low  on  resistance 
(<  200  ohms).  Very  large  MOS  transistors  are  used.  Unfortunately,  the 
large  devices  have  large  parasitic  capacitances.  The  gate-to-drain  and 
.gate-to-source  overlap  capacitances  inject  a  portion  of  the  clock  sig¬ 
nals  into  the  signal  oath.  The  junction  capacitances  ( source-to-bulk  and 
drain-to-bulk)  are  non-linear  functions  of  voltage.  Since  they  shunt  the 
precision  ratioed  capacitors,  they  can  introduce  distortion  into  the 
system.  The  use  of  large  capacitors  can  minimize  both  problems.  An 
upper  limit  on  the  size  of  the  capacitors  i3  set  by  the  clock  period  and 
the  on  resistance  of  the  switches. 
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For  all  the  circuits  except  the  multiplier,  the  feedthrough  of  the 
clock  onto  the  signal  path  results  in  a  D.C.  offset  because  the  same 
amount  of  feedthrough  occurs  on  every  clock  transition.  The  multiplier 
circuit  is  more  complicated  (see  below). 

4..J..  Delay  Line  Limitations 


From  section  3.2,  we  have 


V  .(n+1)  = 
out 


V"> 


If  C^i C^,  the  delay  element  has  introduced  .a  multiol icative  error  in  the 
system  which  must  be  corrected.  This  is  most  conveniently  done  by  the 
microprocessor  which  is  implementing  Durbin's  Recursion  (Apoendix  A). 


Each  delay  element  introduces  a  D.C.  offset  voltage  into  the  delay 
line.  This  is  undesirable  because  the  D.C.  offset  of  each  delay  element 
multiplies  the  signal.  If  the  offset  voltage  is  kept  small  (  <  10  mil¬ 
livolts)  ,  the  contribution  of  this  undesired  signal  feedthrough  will  be 
negligible  due  to  the  effect  of  the  low  pass  filters.  (This  problem  is 
eliminated  if  a  digital  delay  line  is  employed.) 


H-Z'  Multiplier  Limitations 

A/D  Converter  Limitations 

The  A/D  converter  has  3  major  error  sources:  Dloop  gain  error, 
2)loop  offset  error,  and  3)slow  s-’ttline  error. 

Ideally,  the  circuit  has  a  loop  gain  of  exactly  two.  It  is  this 
loop  gain  which  sets  the  radix  of  the  conversion.  The  loop  gain  is  set 
by  the  capacitor  ratios  and  the  open  loop  gain  of  the  op  amps.  On  the 
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breadboard,  an  adjustable  capacitor  was  connected  in  parallel,  with  capa¬ 
citor  "2C"  to  allow  the  loop  gain  to  be  trimmed.  For  the  ID-bit  con¬ 
verter,  the  error  in  the  loop  tain  must  be  less  than  0.12. 


Loop  offset  error  is  caused  by  op  amp  input  offset  voltage,  com¬ 
parator  offset  voltage,  and  clock  feedthrough.  To  insure  less  than  * 
LSB  error  at  the  origin,  we  require 


2  ( V 

loop  offset 

or  V 

loop  offset 


)  <  10.0V 
<  lOmV. 


Slow  settling  is 

caused  by  the 

op 

amp 

outputs 

not 

settling  to 

within  0.12  of  their 

final  value  in 

the 

time 

allowed . 

The 

time  required 

to  settle  has  two  components : RC  time  constants  and  op  amp  settling  time. 
For  this  system,  the  clocks  ran  slow  enough  so  that  adequate  time  was 
allowed  for  settling. 


4.  .2.. 2,.  J3/A  Converter  Limitations 

The  sources  of  error  in  the  D/A  are  Dc.ioaeitor  ratio  error,  2)op 
amp  offset  voltage,  3) finite  op  amp  gain,  4)clock  feed  through,  and 
5) slow  settling  error. 

Asr in  the  A/D  converter,  the  radix  of  the  conversion  is  set  by 
capacitor  ratios: 

radix  s  1  +  7  . 

To  guarantee  *  "J  LSB  linearity,  the  capacitors  must  match  to  better  than 
9.12. 

Capacitor  C£  is  always  discharged  on  £.  If  the  bit  of  the  digital 
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word  is  a  "1",  the  capacitor  is  discharged  to  the  virtual  ground  at  the 
op  amp’s  inverting  input.  If  the  bit  is  a  ”3",  the  caoacitor  is 
discharged  to  ground.  The  op  amp  virtual  ground  is  at  a  potential 
^offset  above  ground.  To  assure  proper  division  of  V^n2  by  factors  of  2, 
C-3  must  be  discharged  into  the  same  potential  on  a  "1"  and  on  a  "0".  The 
op  amp's  offset  voltage  was  externally  nulled  to  ground  potential  to 
eliminate  this  error. 


Due  to  the  finite  gain  of  the  op  amp,  there  is  a  potential  differ¬ 
ence  between  the  inverting  and  non-inverting  terminals: 


V  -V 


out 


open  loop  gain* 


This  difference  gets  larger  as  VQut  increases.  When  the  D/A  converter  is 
processing  the  LSB,  this  difference  (worst  case)  is 


y  _y  -  _ _ 

-  open  loop  vain* 

To  maintain  ±  ~  LSB  accuracy,  the  op  amp  must  have  an  open  loop  gain  of 
at  least  203 0.  For  the  LF356,  the  open  loop  gain  is  about  200,000.  This 
gives  a  difference  voltage  of  only  50  microvolts,  which  is  negligible. 


The  dominant  limitation  in  the  D/A  performance  is  the  clock 
feedthrough.  The  logic  signals  inject  charge  onto  capacitor  C2.  Unfor¬ 
tunately,  two  of  the  logic  signals  are  the  digital  word  and  its  comple¬ 
ment.  Since  the  digital  word  is  always  changing,  the  amount  of 
feedthrough  onto  C2  depends  on  the  digital  word.  This  input-dependent 
feedthrough  is  added  to  the  D/A  analog  output  voltage.  This  error  was 
larger  than  all  other  errors,  setting  the  lower  limit  on  the  analog  out¬ 
put  voltage  at  aoproximately  5  millivolts. 


To  make  the  D/A  converter  compatible  with  the  McCharles*  A/D  con¬ 
verter,  it  ran  at  the  same  clock  frequency  (95  kHz).  Therefore,  the  D/A 
converter  had  adequate  time  for  settling. 

Filters 

Tne  simplicity  of  the  filters  results  in  few  probLem3.  The  op  amp 
input  offset  voltage  and  clock  feedthrough  result  in  a  D.C.  offset  vol¬ 


tage  at  the  output.  The  effect  of  non-linearities  due  to  junction  capa¬ 
citances  in  parallel  with  the  switched  capacitors  were  minimized  by 
using  large  capacitors. 
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5..  Measured  Results 

The  multiplier  was  tested  as  a  unit.  First  the  A/D  and  D/A  convert¬ 
ers  were  checked  for  proper  operation,  and  then  the  multiolier  was 
tested.  Figure  13  shows  the  multiplier  operating  in  the  squaring  mode. 

The  multiplier  output  voltage  has  a  range  from  10.07  down  to  5mV 
(see  section  14.2.2).  This  is  a  6odb  range. 

The  output  filters  were  checked  for  dynamic  range  and  frequency 
response.  The  maximum  input  volt3ge  to  the  filters  is  10.0V.  The  noise 
level  in  the  filters  is  below  1  millivolt.  This  corresponds  to  a 
dynamic' range  of  better  than  80db. 

The  frequency  response  for  output  filters  k=0  and  k=9  are  shown  in 
figure  1U.  Notice  the  similarity  in  the  characteristics.  The  zero  in 
the  transfer  function  seems  to  have  little  effect. 

A  breadboard  system  similar  to  the  one  described  here  has  been  used 
to  process  speecn  t°1.  Figure  15  shows  the  result  of  fitting  a  q-pole 
L.P.C.  transfer  function  to  the  sound  / i/  (the  ’e’  sound  in  ’beet’).  The 


fit  is  very  good. 


3.-  £ uturc  Conslderat ions  one!  duv 'rested  I  morovements 

The  system  proposed  here  is  feasible,  but  is  not  the  best  aDoroach 
to  oe  taken  for  an  integrated  version.  An  integrated  version  could  bene¬ 
fit  from  a  number  of  changes: 

[1]  A  digital  delay  line  should  be  used  rather  than  an  analog  delay 
line.  This  is  more  accurate  and  requires  less  silicon  area. 

[2]  One  MDAC  can  do  the  work  of  ten.  The  one  MDAC  could  be  time  multi¬ 
plexed  to  handle  all  10  multiplications  in  the  125  microsecond 
clock  period  (fig.  16).  This  is  not  a  severe  requirement,  and  it 
would  take  less  area.  If  only  one  MDAC  is  used,  then  all  10  auto¬ 
correlation  values  will  have  the  same  scale  factor  error  due  to  the 
gain  error  of  the  A/D  and  MDAC.  But  since  the  solution  of  Durbin’s 
Recursion  depends  only  on  the  relative  values  rathe"  than  the  abso¬ 
lute  values  of  the  autocorrelation  function,  this  error  need  not  be 
corrected . 

[?]  The  bufter  amplifiers  in  the  output  filters  can  be  time  multi¬ 
plexed.  clash  filter  needs  buffers  between  staves.  To  save  area, 
three  buffer  amplifiers  can  be  shared  among  all  10  filters  if  the 
filters  are  driven  by  10  skewed  clocks  (all  at  3kHz). 

[4]  An  A/D  converter  is  required  at  the  output  of  the  system  to  convert 
the  analog  output  voltages  into  digital  words  for  use  by  the 
microprocessor.  If  all  the  autocorrelation  outDut  voltages  are 
free  of  D.C.  offsets,  then  we  can  use  R(0)  as  the  reference  voltage 
of  the  A/D  converter.  (Recall  that  R(0)  J>  RW  for  all  T. )  This 
results  in  an  automatic  scaling  and  normalization  of  tne 
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autocorrelation  values.  This  should  ease  the  work  of  the  micropro¬ 
cessor  . 

[5l  To  allow  the  system  to  handle  a  wide  dynamic  rann;e  of  speech  sig¬ 
nals,  the  L.P.C.  system  should  be  preceded  by  some  type  of 
automatic  ??ain  control  which  would  take  full  advantage  of  the  A/D 
converter.  The  <?ain  factor  would  have  to  be  transmitted  with  the 
other  model  parameters  to  assure  accurate  synthesis. 
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APPENDIX  A 

An  autocorrelation  L.P.C.  is  based  on  minimization  of  the  square  of 
tne  predictor  error  over  all  speech  samples.  If  we  call  e(i)  the  predic¬ 
tor  error  for  sample  "1"  and  s(l)  the  predicted  speech  sample,  then  the 
error  can  be  written  as 

e( i)  =  s( i)  -  s( i) . 

For  an  all-pole  model  of  the  vocal  tract,  the  predicted  speech  sample  is 
?iven  as  a  linear  weighted  sum  of  the  past  p  speech  samples  (p  is  the 
order  of  the  predictor).  In  fact,  the  model  parameters  are  the  wei^ht- 
in?  coefficients: 

P 

s(  i)  =  2  ’.s  (  i-k)  . 

k=1 

The  a,K’s  are  the  model  parameters  discussed  in  section  2  (also  called 
the  L.P.C.  coefficients). 


where  R(k)  is  the  kth  autocorrelation  la:?  of  the  sampled  speech  (see 
section  2,  equation  f?)).  If  we  have  the  p+1  values  of  the  autocorrela¬ 
tion  function,  the  model  parameters  can  be  obtained  by  solving  this 
matrix  equation. 


Fortunately,  the  square  matrix  has  a  special  form  (it  is  Tosplitz). 
Therefore  the  a^'s  can  be  solved  for  recursively.  The  mo3t  efficient 
method  is  Durbin's  Recursion  [3]: 


1. 

2. 

3. 

M. 

5. 

6. 
7. 
3. 


a 


1  _  ALU 

1  '  R(0) 


Ed)  =  [1  -  (aJ)2]R(0) 


for  i=2  to  p 

i*1  i-1 

R(i>  -  Ja  R(i-J) 

ai _ kjJ _ 

i  '  E(i-i) 


1 

V 


i-1 

i-j 


1  <  j  <  i-1 


fi(  i)  =  (1  -  (aJ)2]E(i-1) 
next  i 


end 


LU 

where  refers  to  the  j6n  model  parameter  of  an  i-pole  model.  When  i=D, 
we  have  calculated  the  model  parameters  for  the  p-pole  model. 

Durbin’s  Recursion  requires  about  100  multiplications  (or  divi¬ 
sions)  and  about  MOO  additions  (or  subtractions)  for  a  9  polo  model. 
With  the  faster  microprocessors  (i.e.  Texas  Instruments  9900),  the 


recursion  can  be  completed  with  the  transmission  interval  (t  10  mil¬ 
liseconds).  Also,  Durbin's  Recursion  uses  only  18  storage  locations. 
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Introduction 


The  design  of  integrated  monolithic  frequency  filters  has  made  rapid 
advances  in  the  last  few  years.  Analog  sampled  data  responses,  which 
are  continuous  in  amplitude  and  discrete  in  time,  have  been  successfully 
implemented  in  NMOS  switched  capacitor  technology,  and  numerous  fixed 
coefficient  filters  have  been  fabricated  [1]— [4] . 

A  simplified  diagram  of  a  speech  synthesis  system  is  shown  in 
Figure  1.  The  focus  of  this  paper  is  on  the  time  varying  filter  which 
models  the  vocal  tract,  the  response  of  which  must  change  to  match  speech 
characteristics.  The  switched  capacitor  technique  lends  itself  to  pro¬ 
grammable  filters,  since  the  capacitors  which  control  the  vital  charac¬ 
teristics  of  the  response  can  be  made  variable  arrays  by  switching  in 
different  amounts  of  capacitance  rather  than  being  of  fixed  size.  Using 
this  technique  an  electrically  programmable  NMOS  second  order  switched 
capacitor  filter  has  recently  been  developed  by  D.  J.  Allstot  [5].  The 
vocal  tract  can  be  characterized  by  a  ten  pole  response  [6],  [7],  and  a 
speech  synthesizer  can  be  implemented  by  cascading  five  of  these  second 
order  sections.  An  alternative  is  a  single  ten  pole  filter.  Several 
approaches  for  synthesizing  the  tenth  order  transfer  function  with  a 
single  switched  capacitor  circuit  are  discussed  in  this  paper. 
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Chapter  1  -  The  General  Problem 

Before  proceeding  to  specific  synthesis  techniques,  a  general  over¬ 
view  of  the  problem  and  possible  approaches  to  solution  will  be  reviewed. 

The  z-domain  sampled  data  transfer  function  of  the  vocal  tract  can 

10  -i 

be  given  by  T(z)  -  i/(i+£  h  z  ) .  The  h . 1 s  are  generally  called  the 

i=l  1  1 

predictor  coefficients  and  will  be  assumed  to  have  been  obtained  from  a 
previous  linear  prediction  analysis  of  the  speech  [6],  [7].  A  typical 
sample  of  speech  is  used  in  this  paper  to  test  the  techniques  developed. 
Several  basic  filter  structures  will  be  investigated  and  compared  with 
respect  to  sensitivity  of  the  response  to  the  number  of  bits  accuracy 
used  in  realizing  the  filter  coefficients.  Bits  are  used  for  comparison 
because  in  an  IC  implementation  the  filter  parameters  are  implemented  as 
capacitor  ratios  which  are  varied  by  connecting  or  disconnecting  elements 
in  capacitor  arrays.  Binary  weighting  of  the  arrays  is  one  of  the 
easiest  schemes  to  implement.  A  different  weighting  may  be  more  effi¬ 
cient  in  certain  cases,  but  a  binary  code  gives  general  sensitivity  pro¬ 
perties  and  also  a  direct  indication  of  chip  area  required,  since 
increased  accuracy  implies  larger  capacitance  area.  The  upper  limit  on 
accuracies  achieved  is  about  10  bits  [8]. 

The  original  processing  of  the  speech  to  obtain  the  predictor  coeffi¬ 
cients  was  done  at  an  8  KHz  sampling  rate,  so  the  sampled  data  filter 
simulations  will  also  be  done  at  the  same  rate,  with  the  exception  of  a 

few  cases  in  Chapter  4.  To  obtain  the  frequency  response  to  be  duplicated, 

10  -i 

consider  the  given  transfer  function  T(z)  =  1/(1+  2  ^.z  ).  In  this 

i=l  1 

paper  only  magnitude  response  is  considered  since  the  ear  is  not  particu¬ 
larly  sensitive  to  phase.  To  determine  the  frequency  response  from  the 
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z-transform,  T(z)jz  -  exp(jwT  )  is  used,  with  T  =  1/8  KHz. 

s  s 

10 

T(exp(jwT  ))  =  1/(1+  £  h  exp(-jiwT  ))  (1) 

S  i*l  1  s 

|T(exp(jMTa))|  -  1/Zt^  +  (2) 

TR£Al  “  1  +  cos  <oTg  +  cos  2toTg  +  . . .  +  h^g  cos  10wTg  (3) 

"^IMAG  ^1  S^n  “Ts  ^2  S^n  2wTs  +  ...  +  hlQ  sin  10a)Tg  (4) 

|T(exp(juTs))|  -  -  20  log  +  T^J  In  dB  (5) 


The  program  for  calculating  and  plotting  the  magnitude  frequency 
response  at  100  points  between  100  Hz  and  4  KHz  is  given  in  Figure  2. 

The  predictor  coefficients  are  read  into  the  array  H(10),  and  Y(100)  is 
the  output  magnitude  plotted  vs.  the  frequency  array  WW(100)  in  subrou¬ 
tine  PLOT.  The  user  can  specify  the  vertical  range  of  the  plot  by 
supplying  YLIT  and  YBIG  and  setting  IPLT  *  0,  but  if  IPLT  3  1  as  in  the 
example  at  hand,  the  subroutine  calculates  YLIT  and  YBIG  from  the  array 
to  be  plotted,  Y.  The  computer  output,  which  is  the  magnitude  response 
to  be  the  standard  for  all  filters  in  this  paper,  is  given  in  Figure  3. 
From  the  plot  four  distinct  resonant  peaks  are  identified,  and  the  center 
frequency  fg  and  Q  of  each  are  as  follows: 


607  Hz 

qi  * 

5.2 

1543  Hz 

Q2- 

29.7 

2401  Hz 

Q3  = 

8.4 
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3454  Hz 


Q4  -  19.0 


As  will  be  seen  in  a  later  section,  the  thrid  resonance  is  created  by 
two  adjacent  complex  pole  pairs,  but  only  one  center  frequency  and  Q  can 
be  identified  from  the  plot.  The  above  numbers  will  be  the  standards 
of  comparison  in  sensitivity  studies  of  the  various  filters  to  be  inves¬ 
tigated.  On  the  plot  it  is  also  important  to  note  the  heights  of  all 
the  peaks,  as  several  filters  in  chapter  4  realize  nearly  correct  center 
frequencies  and  Q's  but  have  incorrect  amplitude  levels. 

Sampled  data  filter  circuit  elements  (delays,  multipliers,  and 
adders)  will  be  used  for  synthesis  and  simulation.  Typical  examples  of 
switched  capacitor  versions  of  these  circuit  elements  are  given  in 
Figure  4.  Further  discussions  can  be  found  in  [l]-[4],  [9]. 
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2  i,.t=lJ+245 

CALL  PLOT'  LJU.V.  OPTS.  IPLT.  VL  IT.  YB  IG  - 
El  ID 


SUBROUTINE 

PLOT' 

: I.  Y.NPTS.  IPLT.  iT.IT.  BIG' 

I' I MENS  ION 

..  1  500 .' 

.  ,’'  500;  .NA'  26  '  -  i  STC'i  '  , 

DATA  kstp 

4M  ♦ 

.  4H  t.  ,4H  K  .41-  * .  4H 

IF'  IPLT.EO 

.  0  1  GO 

TO  300 

YL IT= 1' 1 1 
LO  100  1=2 

. :  IF  TS 

020 

:L :t=hmini 

'  1  L  .  1  # 

*  \  I  l  1  1 

BIG =Y' 1  1 
10  200  I  =2 

1  ip  .  - 

200 

B  I G  =hMA:  1 1 

'  YE:IG. 

.'I" 

300 

Y  nf  I»jE  "  iC  1  1  n 

-.1  1  r 

do  101  ;--i 

•  l  i 

10 1 

<  I  1  =YL  IT 

+ .  I-  r 

f  F;h  f  1  b  E .  1 

440 

FORMAT  <  IH8 

.  1 4; ; .  1  i  f  1 0 . 3 .1 

Pr'  [NT*  4-iO  . 
PRINT  480 

r  1 

480 

FORMAT'  20; 

.  2H-  I . 

10  t  1  OH - -  :  .  IN-) 

DO  tOO  I = I 

.  2  h 

coo 

MAC  1 1  =I'STP 

«;5  1 

DO  700  i  =  I 

.  NF'T’S 

J-K 

IU=' V'(  Ji-VLIT'  PhIIGE*  100.  +1.5 
IF  'i  ILL  GT.  10  l  1  UJ=  10  I 
IF' IU.LT. H  IU= 1 
ITEM*'  I  Li-  1  '  4 

Mh  '  I TEI M  1  '  =1  STP  IU-4  »■  I TEI T' 

PRINT  490 

490  FORMAT '  '  '  . 201 1.  1H.  .  10 <  10H  .  '  ■ 

PRINT  040 . .  >  J  '  .  •  IIP '  1 .1 .  I  =  1 . 2b> 
b40  FORMAT  *  ‘.4  .F12.4.3M.  1HI.25h4.h1.  1HI) 
TOO  Nh  '  I  TEI  1+ 1 ''  =1  STP '  5"' 

PRINT  48G 

■50  FORMAT.  lh'2.  14'  .  1  1 F  1 0 . 7  > 

PRINT  450. W 
RE TURN 
END 


p'l^u.rg,  7-  *  Frequency 


4-r  C^cu‘i  Elements 


Chapter  2  -  The  Singly  Terminated  Ladder 
One  of  the  simplest  structures  from  a  synthesis  standpoint  is  the 
singly  terminated  ladder  described  by  Mitra  and  Sherwood  in  "Digital 
Ladder  Networks,"  [10].  In  this  paper  a  transfer  function  in  positive 
powers  of  z  is  required,  which  is  obtained  by  multiplying  the  numerator 
and  denominator  of  T(z)  by  z^  so  that 

T(z)  -  z10  x  (l/(z10  +  E  h  z'i+10))  (6) 

i-1 

i  10  i 

since  |z  |  *  1,  the  problem  is  reduced  to  synthesizing  the  equivalent 
transfer  function  (with  respect  to  magnitude  response)  of 

10  q 

T(z)  =  l/(zAU+h1z',+  ...  +hgz+hlQ)  (7) 


The  method  of  Mitra  and  Sherwood  is  directly  applicable  to  the  case 
at  hand,  and  through  a  proper  change  of  variables  the  pure  delays  of 
their  procedure  can  be  replaced  by  Direct  Discrete  Integrator  (DDI) 
sections  shown  in  Figure  4(d). 

As  described  in  the  Mitra  and  Sherwood  paper,  the  general  form  of 
the  filter  without  terminations  is  that  of  Figure  5.  If  the  network  is 
considered  a  2-port,  then  the  terminal  characteristics  can  be  related  by 
a  "C"  matrix  as  follows: 


If  the  ladder  is  terminated  at  the  output  by  letting  A2  ■  -  B^,  then  the 
"C"  matrix  takes  the  following  modified  form: 
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C11  C12 


C21  C22j  rB2 


Expanding  the  second  equation. 


B2  "  C21A1-C22B2’  °r  B2^A1  =  G21^1"K:22^ 


For  a  continuous  time  LC  ladder  network,  the  terminal  voltages  and 
currents  can  be  related  by  the  small  signal  short  circuit  admittance 
parameters  as  follows  [11]: 


i-lCs)  yH(s)  yl2(a)  vl(s) 


i2(s) 


y21(s)  y22(s)  v2(s) 


(11a) 


If  the  network  is  singly  terminated  with  one  ohm  at  the  output,  then  the 
small  signal  voltage  gain  is  given  by 


Av(s)  *  -  y12(s)/(l+y22(s))  (lib) 

In  the  continuous  time  case  the  Cauer  synthesis  technique  is  as  follows: 
Given  Av(s)  *  1/D(s): 

(1)  for  the  even  order  case, 


y22(s)  -  Even(D(s))/Odd(D(s)) 


or  l/y22(s)  =  1/ (Even(D(s))/Odd (D(s) ))  (1 

(2)  Perform  a  Cauer  I  expansion  (for  the  even  order  case)  of 
l/y22(s)  to  obtain  the  element  values  of  the  LC  ladder.  The  resulting 
Cauer  I  expansion  is  of  the  form 
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l/y,,(s)  = 


<io8+' 


V  + 


CgS  +  . . .  + 


C„s  +  ; - 

2  L^s 


(3)  The  desired  A^Cs)  is  realized  up  Co  a  multiplicative  constant 

which  can  be  determined  by  comparing  the  original  Av(s)  and  the  transfer 

function  A^(s)  actually  realized  at  any  convenient  frequency.  DC  is 

often  the  easiest.  Then  a  multiplier  G  »  A  (0)/A'(0)  is  attached  at  the 

v  v 

output  to  complete  the  realization. 

Equation  (10)  suggests  that  a  similar  procedure  can  be  used  for 
calculating  the  coefficients  of  the  sampled-data  ladder  if  is  the 
output  voltage  and  is  the  input  voltage.  A  Cauer  I  expansion  of 
C22  *  (B2^A2^IA1  “  0  can  be  put  in  the  followlng  form: 


a^z  + 


c^z  +  . . .  + 


aQz  + - 

9  a10z 


By  longhand  analysis  of  the  filter  in  Figure  5, 


T£+  . . .  + 


T  +  — = 
9  T, 


Comparing  equations  (14)  and  (15),  the  sampled  data  ladder  requires 


1/T^  -  z~X/a^,  i  -  1  through  10.  Then  the  singly  terminated  ladder 
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takes  the  form  in  Figure  6. 


The  filter  coefficients  are  given  by  P^=  -l/a^,  P2  =  1/a^,  .... 

Pg  *  -l/c^.  P^  *  1/a^,  where  through  ci^q  are  the  results  of  the 

10  9 

Cauer  I  expansion  of  Cj2  “  Even(D(z)) /Odd(D(z)) .  D(z)  *  z  +  h^z  + 

...  +  hgZ  +  h^  as  given  previously.  The  gains  of  the  desired  T(z)  and 

the  sampled  data  ladder  at  DC(z=l).  G  =  A  (z=l)/A  (z=l), 

desired  vrealized 

1° 

where  A  (z»l)  *  T(z)jz  =  1  =  1/(1+  £  h  )  and  A  (z=l)  is 

desired  i=l  Vrealized 

obtained  by  a  longhand  analysis  of  the  ladder  at  z  =  1.  It  is  found 
that 


(z« 1) 


realized 


F1F2F3F4F5F6F7F8F9F10 


(16) 


with 


F1  =  1/(1+^),  F2  -  l/CFj+a^,...,  Fw  =  l/(F9+a1Q) 


(17) 


The  program  used  to  calculate  G  and  P^  through  P^q  is  given  in 
Figure  7.  In  the  main  program,  the  f^'s  are  read  in  to  the  array  H(10) 
and  arrays  A  and  C  are  filled  with  even  and  odd  elements  of  H,  respec¬ 
tively,  plus  A(6)  *  1  to  account  for  the  unity  coefficient  of  z^  in 
D(z).  Then  the  subroutine  COEFFS  is  used  to  calculate  the  Cauer  expansion 
of  Even (D (z )) /Odd (D(z)),  and  returns  the  coefficients  in  array  ALPH. 
Finally  the  filter  parameters  P^  through  P^q  are  computed  as  discussed 
earlier  and  put  in  the  array  PARAM(IO),  and  the  gain  factor  GAIN  is  cal¬ 
culated  in  the  subroutine  GAINP.  The  algorithm  used  in  COEFFS  is 

similar  to  the  Routh  tabulation  in  stability  analysis  of  control  systems 
112]  and  is  illustrated  in  Figure  8,  with  equations  in  the  subroutine. 

The  first  two  rows  of  the  chart  are  the  coefficients  of  the  even  and  odd 


powers  of  z  in  D(z),  respectively,  in  descending  order. 

The  output  of  the  program  for  the  specific  case  given  in  Chapter 
1  is  shown  in  Figure  9.  To  test  the  synthesis  algorithm,  the  Purdue 
digital  simulation  program  DINAP  was  used  [13].  The  frequency  response 
from  100  Hz  to  4KHz  is  shown  in  Figure  10.  Excellent  agreement  with 
the  desired  response  is  observed.  The  results  of  the  sensitivity  study 
made  by  changing  the  number  of  bits  accuracy  in  the  filter  coefficients 
are  shown  in  Figure  11.  For  the  general  sensitivity  analysis  all  para¬ 
meters  were  changed  simultaneously;  no  attempt  was  made  to  vary  only 
selected  multipliers.  To  round  a  number  p  to  n  bits,  p  x  2n  was  rounded 
to  the  nearest  integer  and  divided  by  2n. 

p|n  bits  *  Integer (px2n)/2n  (18) 

The  number  of  bits  required  depends  on  the  accuracy  desired  in  the 
response.  For  a  requirement  of  no  more  than  10%  error  on  all  center 
frequencies  and  Q's,  11  bits  of  accuracy  were  needed.  Below  11  bits 
the  response  varied  randomly  and  for  other  speech  responses  than  the  one 
simulated  here,  the  trends  would  be  different  but  unpredictable. 

The  high  bit  requirement  severely  limits  the  usefulness  of  the  singly 
terminated  ladder  due  to  the  large  chip  area  required  in  fabrication  plus 
the  excessive  amount  of  control  circuitry  needed  to  program  the  filter 
coefficients.  However,  the  concepts  developed  will  be  useful  in  other 
algorithms. 

Despite  the  sensitivity  limitations  of  this  circuit,  it  is  of 
general  interest  to  synthesize  the  filter  with  the  DDI  sections  in 
Figure  4(d).  An  alternative  would  be  a  structure  using  Lossless 


Discrete  Integrator  (LDI)  sections,  which  has  a  transfer  function 


-1/2  -1 

T(z)  ”  z  / (1-z  )  and  differs  only  in  switch  phasing  from  the  DDI 

section  in  switched  capacitor  implementation  [4].  However,  the  synthesis 
procedure  for  the  LDI  transformation  is  theoretically  much  more  diffi¬ 
cult  and  has  not  yet  been  developed. 

The  general  DDI  section  without  subtraction  and  multiplication 

included  is  given  in  Figure  12.  The  transfer  function  V  (z)/VJ  (z) 

out  In 

=  z  1/(1-Z  *■)  is  defined  as  a  new  variable  w  \  so  that  the  previous 
synthesis  technique  for  ladders  with  pure  delays  can  be  used. 

z  ^/(1-z  ^)  *  u  \  or  z  =»  ui  +  1  (19) 


If  z  is  replaced  by  u>  +  1  in  T(z),  all  expressions  are  multiplied  out, 

and  the  coeffieicnts  of  each  power  of  ui  are  collected,  then  the  original 

synthesis  procedure  can  be  used  for  T(u>)  *  T(z)|z  **  u  +  1,  the  transfer 

function  of  the  DDI  filter.  No  simple  relation  like  Equation  (19)  can 

-1/2 

be  derived  for  the  LDI  section.  In  that  case  V  (z)/V  (z)  =  z  / 

out  In 

-1  1/2  -1/2 

(1-z  )  and  u)  *  z  -  z  ,  which  greatly  complicates  the  above  sub¬ 

stitution  and  expansion  process.  The  general  form  of  the  DDI  structure 
is  shown  in  Figure  13.  A  different  DC  gain  factor  needs  to  be  calculated 
also. 

To  find  the  coefficients  of  each  power  of  the  new  variable  ui,  the 
following  procedure  is  used: 
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(w+1)  +h^(urbl)  +  ...  +  hg(crH)+h^ 


*10  9 

u  +gjW  +  •  •  •  + 
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By  long  multiplication  and  collection  of  powers  of  w. 


gx  -  10  +  h1 

g2  “  45  +  9h^  +  h2 

g3  -  120  +  36hL  +  Sh2  +  h3 

g4  ■  210  +  84h^  +  28h2  +  7h3  + 

g5  -  252  +  126hx  +  56h2  +  21h3  +  6h4  +  h5 

g6  -  210  +  126hx  +  70h2  +  35h3  +  lSh4  +  5h$  +  h6 

g7  -  120  +  841^  +  56h2  +  35h3  +  20h4  +  10h5  +  4hfi  +  hy 

gg  -  45  +  361^  +  28h2  +  21h3  +  15h4  +  10h5  +  6h6  +  3h?  +  hg 

g9  -  10  +  91^  +  8h2  +  7h3  +  6h4  +  5h5  +  4hfi  +  3h?  +  2hg  +  hg 

g10  "  1  +  hl  +  h2  +  h3  +  h4  +  h5  +  h6  +  h7  +  h8  +  h9  +  h10  (22) 

The  computer  program  used  to  calculate  the  DDI  filter  coefficients 

and  gain  multiplier  is  given  in  Figure  14.  The  same  subroutine  C0EFFS 

is  used  to  calculate  the  Cauer  I  expansion  and  was  given  in  Figure  7.  The 

arrays  A  and  C  are  modified  with  the  relations  of  Equation  (22),  however. 

The  gain  factor  is  calculated  in  subroutine  GAIND.  The  DC  gain  desired 

10 

is  T(z)|z  ■  1  ■  1/ (1+  2Z  *0  an<1  t*ie  Sain  the  DDI  filter  in  Figure 

i=l  1  10 

13  is  seen  to  be  one,  so  the  gain  factor  is  G  =  1/(1+  ^  h.).  The  output 

i=l  1 

of  the  DDI  calculation  program  is  in  Figure  15. 

The  DINAP  simulation  with  no  rounding  of  coefficients  is  shown  in 
Figure  16.  Again  the  agreement  with  the  ideal  response  is  nearly  perfect. 
The  results  of  the  sensitivity  study  are  in  Figure  17.  Overall  the 
results  are  worse  than  those  of  the  filter  with  simple  delays.  For  10% 
accuracy  on  all  center  frequencies  and  Q's,  13  bits  of  accuracy  are 
required.  Belov  13  bits  the  variations  are  again  random,  and  in  several 


inT*. 
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l»]»] 


DDI  singly  terminated  ladder  is  not  a  feasible  solution  to  the  problem 
either. 
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D  I  MENS  IGH  A  <  6 )  .  C  '  5  ' .  ALFA  (1 Q ) .  PflRAMf  1 0) .  H <  1 
DO  I  1=1.10 

1  READ  2. IN  I ) 

2  FORMAT'  F 5  .  O  ' 

A' 6) =1 

A'. 5)  =H i  2'1 
At4>  >HU) 

A  i.  S'1  «H(  6') 

A<2>  *H  1  8) 

H 1.  1 )  =H  (  10) 

0(5) »H  <  1 ) 

C  '14)  =H  *3) 

L ■ 3 1 «H(5» 

C  ',2)  «H  <7'i 
C  (  1 =H  i.9) 

CALL  C  OEFFS  C  A . C . ALPH ) 

PAPAI-K  1  i  =- l.-'ALPHC10; 

F  ARAM ( 2 )  ■  1  ALPH  < 3 ) 

PhPhM i 3 1 =- 3  - ALPH 1 S ) 

PARAMi.4)  =  l/ALPH':7) 

PhRAMi  5  1  *- 1  'ALPH' A' 

PAF'hCI'  b >  =1- ALPH'- 5) 

P ARAM '  7  >  =-l.  ALPH'  4i 
PAPAM'  3"'  =  1 -  ALPH 1 3 1 
PAP  AM '  9  '  =- 1  ALPH'2‘1 
PAR  AM 1 0  ■'  =1  ALF'H  1  1  ' 

CALL  GA 1  HP  •' ALPH.  H.  GAIN) 

PRINTS 

5  FORMA  I'' ’O' LADDER  PARAMETERS  1  THRU  10'  ) 
DO  3  1  =  1  •  10 

3  PPINT  4.PARAM1  P 

4  FORMAT' '  ' .ElF. 10) 

PRINT  6- GAIN 

b  FORMAT'. ' 0 ' , '  GA  I H  FACTOR  = ' .  E 1 7 .  1 0 ) 

END 


pi  ~1  o  C<^1 


hA  !  F.) .  CC  '"5)  .  HLPHFl 


;•  i !  'EM5  ! 0!  I  0 ' 5  1  •  P  *  4  • .  55  ■  4  ■  .  r  T' :  3 :  •  U 1  3  i .  7V 1  3  • .  i.  •  >'  3 :  . 
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i  T  4  2 '  = '  1  S3  ■  1  :■  t  P '  3  i  '  -  4  P  ■  1  1  *  S3  ( 3  1  •  SS 1  l  < 

T1  •31*1 4  SS 1  1  ’•  4  P  r  4)  >  -  •'  P  ■  1  >  •  SS  <  4  ■  t  55  f  i  ■ 

Ij  ;  i  =  i  i  yT  i.  ;  i  i-s:S  i  2  1  •  -  •  S  S 1  1  1  4  TT 1  2  1  :  ■  TT  ■.  1 ) 

'  J  ■  2  •  *<  '  TTi  1  •  4  55 -3;  •  -  'SS '  1.'  *T7>  Z> )  .•  TT'  1) 


U 1 31 *SS <4 > 

TV 1 )  *  t  (.IJ;.  1 )  »-TT 1  2)i-  C  TT  i  1 ;  +IJ  >.  2  >  •  '•  -  IJ  <  1  ) 

V'-'  4  2 )  *  •  <  IK  1  >  ITT  i.  3 1  i  -  •  TT  4  P  nj  <  5 )  <  >  ■  Ij  <  1 ) 

u>  n  =  4  4 w c  n  4-u 4 2 :•  :> - c u <  p  *w >: 2 >  •< ;< -•  vv d.< 

UK  2)  =U  4  3) 

v* i  4I..J4  P  *W(2.'  1  - 1  VV '  1 >  41.J1  ;  1  '■ '  !,ir  I 

V =UC  2 1 

hLPHA  4  1  j  *Af!  PS)  .  CC  1' 51 
hLdHH  '2;  =CC  4  5  1  Q'T) 

HLPH4H  4  3 )  *Q  4  1'  P'  P 
mlPHh  1  4.'  =R  4  P  SS1  P 
hLPhh (5) =SS (  ! )  TT '  P 
riLPHA 4 5 )  =TT 4  1 IJ 1  P 
hlPHh 47  '  *IJ.  1  '  VV'  1  ' 
mi.PHh t C >  -  VV 4  1 Mi  1 1 
ALPHA  ('j.'  =  I J '  1 1  \ 
hlfhm'  i0i  =; y 
RETURN 
END 

SUBROUT I  ME  Sri  1 1  (P  <  NLP  . ;  IH  .  Sh  1 1  II  P 

I1 1  lit  NS  ION  hLP  '  10.'  .HH'  10' 

p 1  -  1 , i  i +hLP 1  P  1 

F2=  1  •  '.F  1  +hLP'  2  / 1 

F3  =  l  1 F2+hLP  4  3 1 1 

F4=  P  ''F3+HLP1  4)  i 

F5M  '  F-  1+hI.P '  5  '  1 

55  1  1  F5+HLP  'b1'1 

FT  l  'FGi-hLP'71) 

F  0=1  ' F7+HLP 1 C 1 1 
F9*i  ■  FS+nLP' S' i 
F  10=1  1  FS+i'iLP'  10)  ' 

C 1  TG  N  =  F  P  F  2 1 F  3  +  F  4  4  F  5  •  4  F  43 ' :  T  7  f  F  8  4  F  9  » F 1 0 

CCTGM  T  4  1  i i ' H '  l  '•  +HH 4  21  +illl'  3'1  +NH 1  4.*  +HH(5)-‘HH(G)-‘HH(ri  +HHi'S)  4-HH  ( 2)  + 
1 HN '  10> ) 

CM  I1IN=CCTGN  CHTijII 

RETURN 

END 
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LADDER  PARAMETERS  1  THRU  10 
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-. 164120 16 !4E +01 
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- . 2 46 S3  1 1 652E+03 

-.735381 423  IE-01 
. 2D  1 1 142SC6E+0 1 
.28416233 14E-0 1 
. 22 1 734 1S37E+03 
- , 2S00000008E-0 1 

GAIN  FACTOR =  . 33S4094755E+0 1 
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DOT  S£c,+ior\ 


DIMENSION  ft '.6)  .  C (5)  . ALPHC  10)  »H(  10)  .PAPAM(  10) 

DO  1  1=1.10 

1  READ2.HU  ) 

2  FORMAT (F5.0) 

ACg) =1 .0 

A (5)  =45+9-*H ( l )  +H (2) 

A (4) =210+B4*Hi 1)+28*H(2)+7*H(3)+H<4) 

A ( 3 ) =2 1 0+ 1 26*H ( 1 ) +70  fH ( ? ) +35*H <3)+l5*H<4) +5*H ( 5 )  +H  (b ) 

A  ( 2)  =45+3b*H  <  1 )  +28  (2) +21  *H  ( 3)  +1 5*H  (4)  + 10  *H  ( 5)  +S*H  ( 6)  +3*H  (7)+P 

A (  1 )  =  1 +H ( 1 )  +H ( 2 )  +H ( 3 )  +H ( 4)  +H ( 5 )  t-H (6)  +H(7)  +H (3)  +H O)  +H ( 10) 
C(5)=10+H>:iJ 

C < 4) = 1 20+36*H ( 1 ) +P+H C 2 ) +H ( i ) 

C<3) =252+ 1 26+H ( 1 ) +56+H < 2 ) +2 1 *H ( 3 ) +6  41 C 4) +H ( 5) 

C ( 2 )  =  1 20+84*H ( 1 ) +56+H ( 2  >  +35+H ( 3 ) +20*H< 4)  + 1 0*H ( 5 ) +4*H C6>  +H  ( 7 ) 

C  ( 1 )  =  1 0+9*  H  ( 1 )  +8+H  ( 2 ;  +7+H  ( 3 )  +S*H  ( 4)  +5  +4H  ( 5 )  +4+H  (6 )  +3+H  ( 7 )  +2*H  (8 ) 
1H(9) 

CALL  COEFFSO.C.ALPH) 

PARAMC 1) =- l.'ALPH< 10) 

P ARAM (  2  )  =  l /ALPH ( 9 ) 

PARAMC3)  — 1/ALPHC8) 

PARAM(4)»1/ALPH(?> 

P ARAM ( 5 ) = - 1 /ALPH ( 6) 

PARAMO)  =1 /ALPH  (5) 

P ARAM (?)»-! /ALPH ( 4) 

PARAMC8)  =  1  ‘"ALPH  C3) 

PARAMO)  =-l/ALPH<2) 

PARAM  ( 1 0)  =  1  .-ALPH  Cl.) 

CALL  GAIND(H.GAIH) 

PRINTS 

5  FORMATC’O'.- LADDER  PARAMETERS  1  THRU  10') 

DO  3  1=1.10  j 

3  PRINT  4- PAR AM ( I )  } 

4  FORMAT!'  '.E17.10) 

PRINT  6. GAIN 

6  FORMAT C ‘O' . "GA IN  FACTOR  =  ' .E17. 10) 

END 

SUBROUTINE  COEFFSCAA.CC. ALPHA)  G\v»*«  1 

SUBROUTINE  GAIHD(HH.GAINN) 

D I TENS  ION  HH( 10) 

CKTGN=1 .0 

ACTGN  =  1  '(  1+HH(  1)+HH(2)+HH(3)+HH(4)  +HHC5) +HH  C6) +HH  (7) +HHO) +HH  O 
1 HH ( 1 0) ) 

gainn=actgn.-cktgn 

RETURN 

END 
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LADDER  PARAMETERS  1  THRU  19 
-.619 17S0035E+00 
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- . 5 145078412E+00 
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-.81 0S3S2333E+00 
.  1032S06464E +01 
— . 1 369362044E+0 1 
.  1356687387E+0 1 
330339930 1E+01 
. 997200000OE+0 1 
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Chapter  3  -  The  Lattice  Filter 
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Another  filter  structure  with  a  simple  synthesis  procedure  is  the 
lattice  filter  with  reflection  coefficients  as  its  multipliers,  classi¬ 
cally  known  for  its  good  sensitivity  properties.  Using  a  simple  recur¬ 
sive  formula,  ten  reflection  coefficients  r^  through  r  can  be  derived 
from  the  predictor  coefficients  h^  through  h^Q.  There  are  many  forms 
of  lattice  filters,  but  one  of  the  simplest,  requiring  only  ten  multi¬ 
pliers  to  realize  ten  poles,  is  discussed  here  [6],  [7], 

Physically  the  reflection  coefficients  characterize  the  behavior  of 
the  acoustic  traveling  wave  in  the  human  vocal  tract,  but  with  proper 
Interpretation,  the  synthesis  algorithm  to  be  described  can  be  applied 
to  any  tenth-order  all  pole  z-domain  transfer  function  P(z).  If  the 
coefficients  of  each  power  of  z  are  collected  so  that  P(z)  has  the 
form 


P(z) 


10  9 

z  -W  **+ 


+  d9Z+d10 


(23) 


then  d^  through  d^  can  be  interpreted  as  predictor  coefficients  and  the 
corresponding  lattice  filter  realization  can  be  derived. 

The  reflection  coefficients  r±  can  be  calculated  from  the  predictor 
coefficients  by  the  following  procedure: 


h*"1  -  (hJ-h^.jJ/d-r^)  (24) 
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for  j  ■  1  to  i  -  1  while  i  ■  10  decreasing  to  1,  where  h^  is  the  pre¬ 
dictor  coefficient  h^.  The  computer  program  for  the  calculation  is 
given  in  Figure  18.  A(10,10)  is  used  for  reading  in  the  predictor 

coefficients  rather  than  H,  and  A(I,10),  i  =  1  through  10,  are  the 
spaces  used.  The  gain  factor  of  the  filter  is  also  calculated  as  will 
be  discussed  later. 

The  total  ten  pole  lattice  filter  is  realized  by  a  cascade  of  ten 
sections,  each  characterized  by  one  reflection  coefficient.  One  of  the 
simplest  sections,  requiring  only  one  multiplier,  is  represented  in  the 
z-domain  as  shown  in  Figure  19.  For  the  first  and  last  sections  the 
structure  must  be  modified.  The  input  section,  characterized  by  r^Q, 
is  shown  in  Figure  20.  The  bottom  left  connection  is  left  open  since 
there  are  no  stages  to  the  left  of  the  first  stage.  At  the  last  stage 
the  proper  connection  is  as  given  in  Figure  21.  Physical  arguments 
govern  such  an  arrangement  and  are  discussed  in  [6] ,  but  the  same 
termination  applies  to  any  transfer  function  with  the  denominator 
coefficients  interpreted  as  predictor  coefficients. 

The  total  filter  is  shown  in  Figure  22.  By  hand  analysis  it  is 

seen  that  the  DC(z*l)  gain  is  one  (without  the  gain  multiplier  at  output) , 

10  10 

and  the  desired  gain  is  1/(1+^  h  ).  Thus  G  =  1/(1+  ^  h  )  and  is 

i=l  1  i-1 

calculated  in  the  program  of  Figure  18.  The  computer  output  giving  the 
reflection  coefficients  and  gain  factor  is  shown  in  Figure  23.  The 
DINAP  simulation  with  no  parameter  errors  is  in  Figure  24,  and  the 
results  of  the  sensitivity  study  are  given  in  Figure  25.  With  no 
coefficient  errors,  the  agreement  of  the  lattice  filter  response  with 
the  desired  response  is  excellent. 
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The  sensitivity  properties  of  the  lattice  filter  are  much  better 
than  those  of  the  singly  terminated  ladder,  as  only  7  bits  of  accuracy 
are  required  to  assure  center  frequencies  and  Q's  within  10%  of  their 
ideal  values.  Another  important  observation  is  that  the  response  does 
not  degrade  as  severely  as  that  of  the  singly  terminated  ladder  for 
lower  number  of  bits,  and  even  for  4  bits  and  less  the  response  is  still 
reasonable.  The  center  frequencies  are  especially  insensitive,  staying 
constant  down  to  4  bits  accuracy. 

If  a  DDI  implementation  were  desired,  one  might  attempt  to  use  the 
same  transformation  z  =  ai  +  1  as  for  the  singly  terminated  ladder.  The 
new  predictor  coefficients  used  to  find  the  modified  reflection  coeffi¬ 
cients  would  be  calculated  by  collecting  the  coefficients  of  each  power 
of  <u  as  done  for  the  DDI  ladder.  In  Figure  26  the  array  elements  B(I,10), 
i  *  1  through  10,  contain  the  modified  predictor  coefficients,  and  Figure 
27  shows  the  new  reflection  coefficients.  To  implement  the  DDI  lattice 
filter,  each  delay  in  the  original  lattice  is  replaced  by  a  DDI  section. 

The  DINAP  simulation  with  no  coeffient  errors  is  shown  in  Figure  28. 
The  response  is  not  as  expected,  and  although  the  four  resonances  are 
still  present,  they  are  not  sharp  and  the  Q's  are  severely  degraded.  It 
appears  that  the  use  of  DDI  sections  causes  a  multiple  zero  at  DC,  and 
the  DDI  lattice  filter  is  not  acceptable. 
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DIMENSION  A( 10. 10) »R( 10) 

DO  1  1  =  1. 10 

1  READ2.AU.  10) 

2  FORMAT (F5 . 0) 

DO  3  1*1.  10 
k  - 1 1  - 1 

R (K) «A<K,K> 

IF'IK.EQ.  1  .* GO  TO  5 
DO  4  J  =  1,1-1 

4  h(J.K-I)  ='.A-  J.IO-A(K.K)^<K-J.K))/Cl-R  00**2) 

3  CONTINUE 

5  CONTINUE 
PRINTS 

3  FORMAT*: ’O' ,  'FILTER  PARAMETERS  1  THPIJ  10') 

DO  b  1=1,10 
S  PPINT7.PU  1-D 
7  FORMAT ( ’  ’.E 17.10) 

GAIN»1/U+AUU0)+A(2.  10)+A(3. 10)+A<4. 10)+A(5.  10)+A(S.  10)+A(7, 10 
1 +A i. 3 .  10)+A(9.  1 0  >  +A ( 10,  10) ) 

PRINTS. GAIN 

9  FORMAT < ’ O ’ . ' GA I N  FACTOR  =  ' . E 1 7 . 1 0) 

END 


l~l  c^are  1%  ’  £e£l€cH«i' 

ce«.  C«. 
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FILTER  PARAMETERS 
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16690979 66E+00 
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. 1781659174E+00 
.  19689606 18E+00 
. 665070045 1E+O0 
- .20 1643768 1E+00 


1  THRU  10 
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Chapter  4  -  Pole  Transformation  for  DPI  and  LDI  Implementation 

The  last  alternative  to  be  discussed  for  synthesizing  T(z)  is 
transforming  the  z-plane  poles  of  T(z)  to  s-plane  poles  and  constructing 
an  equivalent  continuous  time  filter  which  in  turn  is  implemented  using 
sampled  data  filter  techniques.  The  procedure  and  limitation  of  the 
method  will  be  discussed. 

To  find  the  poles  of  T(z),  the  equation 
10  9 

z  +hjZ  +. . .+hgZ+h10  ■  0  (25) 


must  be  solved  for  its  ten  complex  roots.  Using  equations  derived  in 
Hostica’s  paper  [3],  the  five  corresponding  continuous  time  center 
frequencies  and  Q's  are  calculated  as  follows,  where  R  and  0  are  the 
radius  and  angle  from  the  origin  of  each  z-plane  pole.  Since  the  z-plane 
poles  come  in  complex  conjugate  pairs,  only  one  of  each  pair  is  needed 
for  calculation: 

fQ  -  <fg/2TT)/e2+UnR)2  (26) 

Q  -  (-7rf0)/(fs2nR)  (27) 

The  equivalent  continuous  time  transfer  function  T(s),  valid  up  to 
a  gain  constant,  is  given  by 


T(s)  -  n 


2uf, 


i-i  2  -‘0  2 

sZ+(— — -)s+(27rf  )z 
vi  ui 


After  carrying  out  the  multiplication 

1 


T(s) 


10  9 

Sj^s  +s10s  +...+S2S+s^ 


(28) 


(29) 


Standard  s-domaln  procedures  can  be  used  to  synthesize  the  continuous 
time  filter.  A  singly  terminated  LC  ladder,  shown  in  Figure  29,  is 
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chosen  here  because  of  its  simple  synthesis  algorithm  and  the  presence 
•*0  of  design  techniques  for  switched  capacitor  implementation  of  sampled 

data  ladder  filters.  The  basic  algorithm  for  finding  the  L  and  C  values 
has  already  been  introduced  in  chapter  2  and  is  summarized  as  follows, 
with  the  basic  ladder  network  now  characterized  by  its  open  circuit 
impedance  parameters  [11]: 


r  — 

1 

r 

vL(s) 

su<9> 

z12<9) 

ix(s) 

v2(s) 

_  J 

s2l<s) 

z22(a> 

i2(s) 

(30) 


Now  the  small  signal  voltage  gain  is  given  by 

Av(s)  -  z^sJ/U+z^s))  (31) 

(1)  Divide  the  denominator  of  T(s)  into  even  and  odd  parts 


T(s)  -  l/(Deven(s)+Dodd(s)) 


(32) 


T(s)  *  z10(s)/(l+z11(s))  =  (1/D^jj ) / (1+D  (s) /D  ,,  (s) )  (33) 


(2)  Put  T(s)  in  the  form 
!12(s)/(l+zuUy  -  v-'“odd(s)"  "even""'  ‘'odd' 

(3)  With  a  one  ohm  source  resistance,  z^(s)  is  the  impedance  looking 

to  the  right  of  the  resistor,  and  the  L  and  C  values  are  found  by  a  Cauer 

1  expansion  of  z,-(s)  =  D  (s)/D  ,,(s).  The  form  of  the  expansion  is 
1J.  even  odd 

(35) 


*n(s) 


D  (s)/D  ,,(s)  =  L.s  + 
even  odd  1 


C2S  + 


L^s  +. • •+  ~ 


L.s  +  — 

9  V 


(4)  By  inspection  of  Figure  29,  the  DC  gain  of  the  continuous  time 


ladder  is  one,  so  to  make  the  DC  level  the  same  as  T(s),  a  voltage 

10 

divider  of  value  G  *  l/(l+£  h^)  is  needed  at  the  output,  as  shown  in 


i-1 


Figure  30.  It  is  important  to  note  that  the  gain  adjusting  network 
must  not  load  the  LC  ladder  or  else  the  basic  frequency  response  will 
be  changed.  This  requirement  is  satisfied  in  a  switched  capacitor 
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implementation  if  the  op  amps  are  nearly  ideal. 

The  computer  program  used  to  calculate  the  LC  ladder  parameters  is 
given  in  Figure  31(a),  with  subroutines  CALC  and  COEFFS  in  Figure  31(b) 
and  the  complex  root  finder  PROOT  in  Figure  31(c).  The  constant  1.0 
(with  multiplies  the  z^  term  in  equation  25)  followed  by  the  predictor 
coefficients  h^  are  read  into  the  array  T  (since  H  is  used  in  PROOT) 
and  reversed  in  order  in  array  A  for  input  to  PROOT.  A(l)  is  the  coefficient 
of  the  constant  term  of  the  polynomial  to  be  solved,  while  A(ll)  multiplies 
the  z^®  term.  The  other  parameters  of  PROOT  are  as  follows:  N  is 
the  degree  of  the  polynomial  to  be  solved.  NPLUS2  =  N+2 

U  and  V  are  arrays  containing  the  real  and  imaginary  parts  of  the  solution. 
Dimensions  of  U  and  V  are  N+2,  but  only  elements  N  and  lower  contain  useful 
information. 

H,  B,  and  C  are  arrays  of  dimension  N+2  used  for  temporary  storage. 

CONV  is  the  difference  required  to  each  root  in  successive  iterations  for 
the  root  to  be  considered  final. 

Next  the  relations  of  equations  26  and  27  are  used  to  calculate  the 

five  continuous  time  center  frequencies  and  Q’s.  Subroutine  CALC  multiples 

out  the  denominator  of  T(s)  and  returns  the  coefficients  of  each  power  of 

s  in  the  array  S,  with  S (11)  being  the  coefficient  of  s10  and  s(l)  that 

of  the  constant  term.  Arrays  AI  and  Cl  collect  the  coefficients  of  the 

even  and  odd  powers  of  s,  and  subroutine  COEFFS  is  again  used  to  do  the 

Cauer  I  expansion.  The  L  and  C  values  are  returned  in  the  array  ALPH, 

with  L1  *  ALPH(l)  and  C.-  =  ALPH(10).  In  Figure  31(c),  the  gain 
1  1  10 

adjusting  resistor  R  =  l/(l+23  h  )  is  calculated.  In  the  program  T(2) 

i=l  1 

through  T(ll)  are  h^  through  h^. 

The  computer  output  of  the  program  is  in  Figure  32,  and  the  SPICE 
simulation  is  given  in  Figure  33.  From  Figure  32  it  is  seen  that  two 
adjacent  complex  pole  pairs  form  the  third  visible  resonance.  As  is 
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characteristic  property  of  continuous  time  filters,  the  correct  resonances 
are  realized,  but  the  amplitude  of  the  response  falls  off  much  faster  than 
that  of  the  discrete  time  response  to  be  duplicated.  It  is  clear  that  a 
discrete  time  implementation  of  the  continuous  time  filter  must  be  used 
to  raise  the  levels  of  the  higher  frequency  poles. 

The  first  steps  in  the  synthesis  of  the  DDI  or  LDI  implementation 
of  the  continuous  time  LC  ladder' are  writing  the  describing  equations 
of  the  circuit  and  constructing  the  all-voltage  signal  flow  graph  with 
inverting  integrators  that  represents  the  equations.  These  steps  are 
shown  in  Figures  34  and  35(a),  respectively.  In  35(a)  the  currents  at 
the  lower  nodes  are  converted  to  voltages  by  multiplication  by  a  scaling 
resistor  R,  chosen  to  be  one  so  that  all  other  gains  in  the  circuit  stay 
the  same.  The  DDI  and  LDI  sampled-data  filters  have  exactly  the  same 
structure  as  the  signal  flow  graph,  but  s-domain  inverting  integrators 
are  replaced  by  their  discrete  time  approximations,  as  shown  in  Figure  35(b). 
To  simulate  half  delays  in  DINAP,  full  delays  are  used  with  the  sampling 
rate  doubled.  The  multipliers  are  given  by 

Pi  "  1/(fs(Li  or  Ci)}  (35) 

In  Figure  31(c)  the  calculation  is  done,  with  FCLOCKD  being  the  sampling 

rate  of  the  discrete  time  filter.  It  will  be  seen  that  changing  this 

parameter  has  a  large  effect  on  the  frequency  response.  For  the  discrete 

10 

time  filter,  the  gain  multiplier  is  simply  G  =*  1/(1+^  h.)  and  is 

i=l  i 

attached  after  the  V  .  node  in  Figure  35(a). 

out 

The  output  of  the  program  for  f  =  8kHz  is  shown  in  Figure  32. 

Similar  runs  were  done  for  f  *  10.88  kHz,  IskHz,  and  120kHz.  The 

s 

computer  outputs  are  not  included  here,  but  the  DINAP  responses  will  be 
investigated  in  detail. 
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First  the  implementation  using  DDI  sections  as  replacements  for 
integrators  is  considered.  The  transformation  is 

s  -  (z-l)/T  ,  or  3  -  sT  +  1  (36) 

3  S 

and  the  block  diagram  is  shown  in  figure  36(a). 

If  any  pole  z^  has  a  magnitude  greater  than  one  (lies  outside  the 
unit  circle),  the  circuit  will  be  unstable. 

|*kl  -  /(Re(skrs)+l)2+(Im(skTs))2  (37) 

For  a  complex  pole  pair  with  center  frequency  and  Q  given,  the  s-plane 
poles  are 


s  -  -irf0/Q  ±  j2irf0/l-l/(4<n  (38) 

For  fg  ■  3462  and  Q  -  19.6  as  given  in  Figure  32  (difference  between  3462 
and  19.6  and  observed  f^  =*  3454  and  Q  *  19.0  in  earlier  plots  is  due  to 
plotting  at  only  39Hz  intervals),  sfc  *  -554.9 ± j21745.3,  and  for 
Ta  -  l/8kHz,  skTs  -  -.06936 ± j 2. 718.  |zfc|  -  /(l-.069362)+(2.7l8)2  »  2.87  >  1, 
so  the  DDI  filter  will  be  unstable.  Thus  in  general  the  DDI  transformation 
is  limited  in  usefulness  and  not  a  possibility  here. 

The  same  filter  can  be  implemented  using  the  LSI  transformation,  for 
which  the  basic  equation  is 

s  -  (z1/2-z_1/2)/Tg  (39) 

with  the  block  diagram  shown  in  Figure  36(b).  The  effects  on  the  frequency 
response  are  different  than  in  the  DDI  case  and  are  derived  as  follows  [14]: 
For  frequency  response,  s  *  ju>  and  z  ■  exp(jX),  where  X  is  the  discrete 
time  frequency  of  the  filter 


ju>  -  (1/T  )(exp(jX/2)-exp(-jX/2)) 
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(41) 


u>  -  (2/T  )sin(X/2) 
s 

f  -  (f  /ir)sin(X/2)  (42) 

8 

When  implementing  Che  continuous  time  filter  with  LDI  sections 

continuous  time  frequencies  f  are  transformed  to  discrete  time  frequencies 

X.  A  plot  of  equation  42  is  given  in  Figure  37.  It  is  evident  that  a 

frequency  expansion  occurs,  with  the  distortion  increasing  with  frequency. 

In  discrete  time  filter  theory  X  »  it  corresponds  to  f  ■  f  /2,  and  the 

s 

baseband  signal  exists  for  |xj  <  ir.  The  basic  limitation  of  the  LDI 

transformation  is  that  since  |sln  X/2 |  <  1,  only  continuous  time  frequencies 

less  than  f  /ir  can  be  seen  in  the  sampled  data  frequency  response.  For 
s 

the  case  at  hand,  to  put  the  3462  Hz  center  frequency  within  the  baseband 

region  (neglecting  putting  it  correctly  at  3462  Hz) ,  the  LDI  filter  must 

run  at  a  minimum  sampling  rate  of  f  =  it  x  3462  *  10.88  kHz.  For  rates 

less  than  that  value,  the  last  resonance  will  be  lost,  and  the  higher  the 

sampling  rate,  the  closer  it  will  be  transformed  to  the  desired  discrete 

time  frequency,  X  ■  (3462/8000)  x  2ir  ■  2.72. 

Unfortunately  as  the  sampling  rate  of  the  filter  is  increased, 

the  response  takes  on  more  of  the  qualities  of  the  continuous  time 

response,  with  excessive  amplitude  rolloff  with  frequency.  DINAP  responses 

of  the  LDI  filter  for  sampling  rates  of  8,  10.88,  15  and  120  kHz  are  shown 

in  Figures  38,  39,  40,  and  41.  With  f  *  120  kHz  the  fourth  resonance 

s 

is  visible,  but  it  is  so  severely  reduced  in  amplitude  that  is  is  not 

effective.  For  f  “  10.88  and  15  kHz,  the  fourth  resonance  is  present 
s 

but  not  clearly  identifiable,  and  the  second  and  third  peaks  are  already 
reduced  from  the  desired  levels. 

In  actual  switched  capacitor  Implementation,  the  path  representing 

-1/2 

the  source  resistance  will  contain  a  half  delay;  R  *  z  instead  of  1. 
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r 


L/ 


i 


Simulations  were  made  of  this  effect  with  the  same  four  sampling  rates, 

and  no  difference  in  the  response  was  seen.  That  result  is  somewhat 

-1/2 

surprising  since  for  frequency  response  z  -  cos(wT  IT)  -  jsin(u>T  / 2), 

S  8 

and  the  source  impedance  becomes  complex  for  low  sampling  rates.  For 
example,  with  w  «  (2ir)(2kHz)  and  T  *  1/8  kHz,  z  -  cos(ir/4)-j 
sin(tr/4)  -  Y  “j  T  * 

In  a  different  application  with  lower  frequencies  of  interest,  a 

singly  terminated  LDI  ladder  could  be  useful,  and  prewarping  of  the  S“plane 

poles  could  be  performed  by  the  relation  in  equation  42  to  obtain  a  more 

correct  sampled  data  frequency  response. 

To  see  the  general  properties  of  the  structure,  a  sensitivity  study 

was  made  of  the  three  center  frequencies  and  Q’s  realized  (although 

incorrect)  by  the  LDI  ladder  for  f  ■*  8  kHz.  The  results  of  the  study 

s 

are  shown  in  Figure  42.  With  full  precision  on  coefficients,  the  discrete 
time  center  frequencies  are  in  error  as  predicted  by  the  relation 

>  -  2  sin_1(irf/f  )  (43) 

s 

For  example,  the  1537  Hz  peak  is  expected  to  be  seen  in  the  sampled  data 
response  at  X  *  2sin  ^(*•1537/8000)  *  1.296,  which  for  an  8  kHz  sampling 
rate  corresponds  in  frequency  to  f  =  (1.296/2ir)  *  1650  Hz,  close  to  the 
1660  Hz  observed.  From  Figure  42,  8  bits  of  accuracy  are  required  to  keep 
all  center  frequencies  and  Q’s  within  10%  of  the  values  with  no  errors  in 
the  coefficients,  better  than  the  singly  terminated  ladder  of  chapter  2, 


but  worse  than  the  lattice  filter. 
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DIMENSION  A>  124 .Vf 12) .H( 12) ,B( 12)*C'124  ,T<  1  P  .UU(5) ,  W'5)  - 
! F0<5 ) . G i 5 ) .  U 1  1 2 •  UG  <5 > .  S  i  1 1 4 ,  A I  i 6 4 .  C  I f  5 4  .PAPAMU0) .  ALPHi  10 
PI-3. 14159264 
FCLOO  *6000.0 
h  ■  1 2  i  -0 . 

M»  10 

NFLUS2-12 
DO  1  1*1,11 

1  PEAD2. T  •  1 4 

2  F  OF'I  II  I  T  <  F5 . 0 1 
DO  5  1=1,11 

5  A 1 1 ) « T < 1 2-  I) 

C  FIND  Z  PLANE  POLE'S 

C ALL  FF'OOT i H .  A .  U ,  V.  H ,  B.  C ,  CONV.  NPLUS2) 

DO  10  1=1.5 
UUi  n=U(2t  I  -  1 ) 

10  VVC  I )  =V(2+  I  -  1 4 
PRINT  30 

30  FORMAT i ’0* , ’Z  PLANE  POLES') 

PRINT  33 

33  FORMAT i ’  '.'REAL  PhRT  IMAG  PAPT' ) 

DO  31  1=1, 10 

31  PRINT  32. U< I) . V( I) 

32  FQRMATC  ' , E 1? . 10. 2X,E 1? . 10) 

C  FIND  S  PLANE  POLES 

DO  11  1=1,5 

D I  V*UU  ( 1 )  ✓  SORT  C  UU  (I )  **2+VV  ( I )  **2 ) 

THETA-ACOSiDlV) 

R  -SORT 1  UU  >.  I  >**2+VV<  I)«2) 

FOi  I ;  * '  FCLOCI'  '(2  !:P  I ) )  SORT (THE  rA**2+AL0G (R)  k*2)  ) 

UO  >.  !  I  =2  t-p  I  >FG '  I  :> 

11  O'.  I)  *1  -F'l'FOi  I) )  '  FCL0CK*ALOG(PJ> 

PRINT  60 

60  FORMAT i  '  0 ' ■  ’  ') 

I'O  40  1  =  1,5 

40  PRINT41. 1.FOiD  .  I.Q'  I) 

41  FORMAT**  *  . '  F  0  ( ’ ,  1 1  ,')=',  E  1 7 .  10,'  QC.Ii.')*' 

C  FIND  COEFFS  OF  EACH  POWER  OF  S 

CALL  CALCOJO.Q.S) 

C  COMPUTE  COEFFS  FOR  ELEMENT  FINDING  SUBROUTINE 
h  l  (6) =Si 1 1) 

A I  '.5 )  =0(9)  J 

A 1 1 4)  *S  1.7)  & 

A I  (3)  =S  i.5)  'A 

A  I f  2  »  =S i 3>  ^  *T 

AIil)*Sil)  <7 

C 1(5) -Si  10)  £  £ 

C I  *  4 '  =S  ■  8 ) 

C I (3) =Si6)  $  p 

CH2)=Si4) 

CHD-S(2) 

C  COMPUTE  ELEMENT  VALUES 

CALL  COEFFS (A I , C I , ALPH)  A? 

PIJNIJ . 2  COMPILER.  MAR  73  VERSION  % 


I. Oil) 

.11-') 


QC.Ii. 


■  E  17 .  10) 


SUBROUTINE 


4-y 


A? 


PRINT  13 

13  FORMAT'. '0',' ANALOG  L  AND  C  VALUES  LI  THRU  CIO') 
DO  20  1  =  1.  10 

20  PRINT  14, ALPH ( I) 

14  FORMAT*'  '.E17.1Q) 


picture.  '^(pT)  * 

c,  Oroya*  -£>'  a*oaU| 
LC  \idd" 


[.  I  MEMO  '.OH  Iji'ii-  _. •  .  I'm.  5  I  1  !  .1 
ST3  -IJOO 1 2  '  NO 1  2  |  HJOO <  P  uO  <  !  > 

■i T 2 »UOO  '£■*-*''  +I.IOO >2t  r*2+< MO <  1  >  ♦tJOO >  2) )  /  t 00  <  1  >  <00  (2)3 
STJ  *  i  UOO  i  111*  Ujiji  2  ’  *  *2  J  3  GO '  1  >  +  ( LJOO  *  1  3**2)  wJQOiS)  >  ''OG  (2) 

ST0-UOO1  1  1  +2  MJOQ  ‘  2 )  i *2 
STS  'IJOO '  4 '  00 ' 4 '  +1.100  <  7  '  00  >  3 ) 

ST7-U00<3'  *  *2+WOU>  4;  i  *2+‘ UOO<33  HOC ''4;  )^i'OU< 3.*  *QQ<4j) 

STb  *  *  IJOO  i  3  1  UOO  ■  j  m-2.1  •  00  *  3  3  +  (  1  UOO  1  3 3  *2  3  *U00  <4>)  ' 00  1  4) 

ST5  =UGO  |  3  ■  **-2  *UOQ  i  4 1  4  1 2 

T 1 =wOO 1 5  *  -  00 1 5* 

TO  »W00 1  5  1  *  *2 

.  m  -sTj  ■  j 

X6-ST2+ST7  +ST3*ST8 

:  5  -ST1  +STb+ST3*ST7-‘-ST2  t-STS 

;  :4»ST0+ST5+ST3*STb+ST2  f  ST r  +ST 1  *bT3 

X3  ’  S  T  3  *S  T  5  +3  T2  T  S  T6  +S  T 1  ♦ST?  +STQ  *ST8 

X2-bT2+ST5+ST  1  t-STO+STO  «ST7 

::i=GTH<ST5+ST0+STb 

:  ,0  -STO  *ST5 

SS<1P*1 

SS<  10)  ■X7+T1 

SS  *1 0  1  =X6+T  1  *X? +T0 

SS(83  *X5+T  1  *Xb+T0*X7 

SS i 7 > =X4+T 1 *X5+T0*X£ 

SS'.b)  =X3+r  1  •■X.4+TQKX5 
SS  «  5 3  »X2+T  1  TV3+T0*X4 
SS  c  4)  »x  1 +T 1  *X2+  r0:«X3 
SS  <  3) -X0+T 1 *X1 +TQ KX2 
SS(2)  »T1  *X0+T04-X1 
S S  <  1  >  ■TO'*'.'  0 
RETURN 
END 

SUBPOUT  I  ME  C  OEF  F  S  <  AA  .  C  C .  hLPHh  3 

D  I  MENS  I  ON  0 ' 5 1'  ,  P  <•  4  3  .  SS  <  4  3  .  TT  <  3  '<  .  U  <  3 3  .  W '"  2 )  .  I, J '  2 )  .  Hft  <  6 ) .  CC  ( 5 ) .  ALPH 
1  « )0> 

0  *  P  ■•  *  CC  1  5  *  *AA (5>  AA  <  6  '  rCC  < 4) )  >  C C  ( 5 ) 

0 1  2  >  ■  1. 1  C  C  1  5 )  *  hh  •'  4 1  1-i.hh  6  '•  *  CC  <  3 ) 3  '  CC  <  5  j  ^ 

0.1.3 >  ='  i  C L  1  5 )  1 1 ih <3  >  1  -  hh  1  b  1  ICC 1  2)  j ) ■  CC  (5)  W 

0 1  43  =  < '  C  C  <  5  '  *  HA '  2  3  '  -  ‘  HA '  0  1  *C  C '  1 3  P  ✓  C  C '  5  3  Jf 

O' 5)  -HH'.  1  '  jp  . 

R  t  p  *  1. 1  Q  (  1  i  ♦C  C  '  4  j  C  '  5 ' 1 0  ( 2  >  >  3  /  0  '  1  '  jjr  / 

P  i  2)  - 1  (0 1  l  i  ii'i.i")  i  - 1  CC  1  5  1  M2  >'  3'1  * 3  /0 1  1  1  V 

p  i.  3  *  ■  *  i  0 1  1  1 ")  C  C  '  2  11  — 1  C  C  1  5  1  *  0 1  4 1  1 )  /  0  *  1  1  ^  3* 

P  i.  4  *  =  i  i  Q 1  1  1  +  CC  <13  i  -  ( CC '  5  1  *  0 ' 5  1 ) )  /•  0 '  1 )  j 

SS  ( 1  1  s  1  '•  R '  1‘  '  0 1. 2)'  -  '0 1  i  '  *  I-'  *  2)  1 )  -  'P '  1  j  -j  £ 

SS  <  2  >  ■  i.  (pi  l  >  *  0 1  3 *  -  <  0 1  l  1  i  P '  3  ■'  1 3  rP  <■  l )  ^  4 

SS < 3  >  ■  1  i  Pi  P  tQ i  4)  j  - i  O'  1  i  »R <4j  *  )  •  P '  1 

SS(4)=Q<5>  4* 

T T '  1  i  =  i  <SS'.  P  IP' 2) )  - '  P 1  1 .»  *SS<'2.) ) )  -•  SS<  1  >  V  ^ 

TT <  2 )  ■  i 1  SS  i  P  *P  (  3 )  1  -  1  P >  1  1  )C.S  i'3) ) )  -  SS  <  1 3  x, 

TT < 33  *'  *  SS <  1 3  i-F' i  4.'  3  -  i  p '  P  i SS ( 43  3  3  -  SS  <  1 3  ^ 

IJ  (  P  ■  ( <  TT  ( 1 3  SS  i  2  3  1  - '  SS 1  1  1 '  P TT '  2  3  3  3  . '  TT  ( 1 ) 

IJ  <  2  3  - 1  i.  TT  1. 1 3  *SS  <  3  3  3  - 1  SS  ( 1 )  "KTT  ( 3  3  3  )  xTT  ( 1 3  J  A? 

LU33-SSU3 

W(P«(  (U(  13*TT'  2  '  ' - (  TT'  1  MlJ'2)))  Ijr  1  i 
VV«  2)  =  <  (IJ  •  P  *TT<33  )->  TT'  P'  HJ(33  3  3  U' 'P 
IJ  i.  P  ■*  ( i  W 1  13*  U '  2 1  '  - '  IJ  i  p  I  \'V '  2 '  3  /W  f  1  ) 

Uk  2  3  *U  I  3  * 

:  >(  I  IJ  C 1 3  ■(>  W  <2)3-1  VV I  1)  NJ I  2 )  3  3  /  IJ  I  1 3 
Y»U(23 

hLF'Hh  i  I  '»  -HH  i  6  >  c r  >  S  ) 

ALPHA'  2  '  -LC  '  5  )  O'  1  '  , —  rju  i\  • 

ALPHA '  3  '  =C"  1 3  P'P  pv«\u^€  * 

ALPHA  <4)  -P  P  -  SS '.  1 1  1  ./  “ 

ALPHA >  5 >  -S3 <:  1  '•  TT >.  1 3  4b* 

ALPHA  lb  '  »T  i P  U  *  1  1 

hlpha  <  t  >  —  i  j  *  r*  -■  vv  <  i '  cmc^Utioiv 

ALPHA'  S  '  i  1 3  IJ<  !  * 

ALPHA' 9  1  =i.i' p 

i'iLPHAi.  IN  1  -X-  Y 

RE'HJRN  . . 


Fv^u te  31  ^  • 

^v\lpfo%2t|n€S  L.  C, 


-iL  t-  1  s_  i  c. - 


.  run  -Li.rLicr  _  jtfrs  ryt  :•  I  , 

ra.ixu.  -  120000.0 

DO  12  1*1.10 

12  PhPhID  1  >  -  1  -  ■  Ft. LOCI  D'liLPtU  I  >  ) 

PP  HIT  15 

15  FOPMhT DIGITAL  LADDEP  PARAMETERS  1  THRU  10' j 
f<0  lb  11-10 
!b  hP  I  NT  ir.PnPHlU  I  < 

17  format • ■  ■ .cir. io 1 

C  FIND  hNhLOG  PLi  I  STOP  HttD  L>  1 G I T  NL  MULTIPLIER  TO  ADJUST  DC  GFlIN 
SUM*T t  2.'  +  T>  3:  +T'  4 j  +  T'  5 )  +T'.6)  +Ti  ?)  +Tl  3 >  +T '  9;  +T(  10)  +T ( 1 1 ) 
P*iUM 


FOCT0R*!  ■  1+SIJ1I' 

PPIMT18.P 

13  FOPMhT t *0* . 'PESISTOP  VALUE*' .E17. 10) 

PR  INT19.  f  hCTijF: 

19  FOPMhT ('O’.' DIGITAL  GAIN  FACTOR*' ,E17. 10) 

END 

SUBROUTINE  PPOOT<  tl .  A .  1 1 .  V.  H .  B.  C .  CGNV.  NPLUS2 ) 

Ml EMSIOrl  O'  IIPLUS2  ■ .  U 1  NPLUS2) .  Vf  IIPLUS2)  •  H 1 IIPLUS2)  .6<NPLUS2) 
1C'HPUJS2' 

COHV's-1  .E-25 


!,4 


Nf 


.  100. 

.  b  *  r 


nc*n+i 

DO  1  1*1.  tic 

1  H  C I ) *0 1  I ) 

p*0. 

Q*0. 

R-0. 

I REV* 1 

3  I F  (  H  *,  1 )  '  4 

2  NC-NC-i 
VIIC  '  =0. 

U'  lie  J  =0. 

DO  1002  1*1- 

1902  h  <  i > -m  « :  - 1  • 

GO  Tu  3 

4  IF' IIC- 1  '5 

5  IF  >IIC-2  >7 

6  P  =-H '  1  •  H  ■  2  ' 

GO  TO  50 

7  IF  till. -3  >9.3.9 

3  P=H>  2'  ■!!'“' 

O  -Ht  1  '  H'  3  ' 

GO  TO  70 

9  IF'.hBS'H'IIC-I) 

10  I  REV—  IPEV 
I1-HC/2 

DO  1  1  1  =  1 .11 

ML=NC+1-I 

F*H'.HL) 

H- ML  '  =1 K  I '' 

11  H '.  I )  *F 

I F 1 0 ) 13. 12.  13 

12  P=3 . 

GO  TO  15 

13  P-P  Q 

q=|./q 

15  IF'R) lb. 19. 16 
lb  R-1.XR 

19  E=5.E- 10 

B  f  MC )  =H  i  NC ' 

Ctiic '  =h  'tie ) 

B'.nc+n  =0. 
c  itlC  +  1 )  *G. 
iip*hc-i 

20  DO  49  J= 1.1000 
DO  21  1 1  =  1. HP 

I  *1 1C  - 1 1 

e- 1 • =h ' i iTp iE '  i+i ) 

21  C  1  I  1 *B'  I ' +P iC 1  I +1' 
ir  .  hRO '  B '  l  '  II'  1  '  '  -  E 


-v'V 

*5.  .4 

#v 

<* 


xV 


H 1 1  IC -hBS  i" H  ( 2 )  'H  ( 1 ) ) )  10.  19.  19 


Fyue  3i(c) : 

CoA'VinUa'il®^  c^ 

b*SvC 


On . 50 , 24 


||»'«|  I 

-  " Vr- ■  .;•. 


2  1  IF 1 1 1  2  1  1  •  23 

22  R-P+l 
GO  TO  30 

23  P»P-6<  1  1  •’  ' 1 

30  I«0  37  il-I.nF 
I-ilC-1 1 

6  ■  J i *H ■ I ) -P +B ■  I  + 1  1 -U *B  I  +2 i 
77  C  •  I  1  =  B 1  I  i  -  P  *  C  1  I  +  1  1  -  0  't-i.  I  +2  > 

I F  •  H '  2 '32. 3!.  32 

31  IF  1  hBS  iB  >  2  1  -  H  '  1  1  ■  F'ZJ.33.34 

3 2  I F  1  hB’3  1  B 1  _! 1  H 1  2  1  1  -E  1  3  j  •  i3 .  i4 

33  IF 1  hOS <  EM  1  H'  1'  <-EV0.F0.34 

34  CBhP«C ' 2 i -B' 2) 

D=C  •  1>  t  12-i‘BhP  i  l  1  4) 

I F '  L> 1 36  •  35 .  36 

25  P»F'-2 

Q*Q  4-1.0+ 1 .  > 

GO  TO  49 

36  P*P+'  B'  2  1  h"  •  3.'  -DM  »  +C !  4j  )  -  D 
U -v+  (-0(2  f  «•£ BhP+B <'  1 )  4cf  i 3 M  'Ii 

49  CONTINUE 
E-EMU 

IF  i.  E-COHV • 20.20. 40 
40  CONV-E 
GO  TO  20 

50  NC-NC-1 

van:  >  -o. 

IF >  I REV'S  1 -52.52 
si  uac.'-i./P 
GO  TO  53 

52  uao-p 

53  DC  54  I-l.HC 

54  ml  1  *B <  I  +  I  1 
GO  TO  4 

20  MCMiC-2 

if (  ipcv > ri.ro. 72 

71  QP-l./O 
PP=P -(0  *2 . 0.' 

GO  TO  73 

72  GP-G 

PP-P/^.Q 

73  F  =  >.  PP  >  r+2  UP 
IF  1  F"> 74.  75. 75 

7d  uac+i )  — pp 

IJUIC)  — PP 
VXNC+l)  -SORT ''  -F 
vaiCj  —  van;+n 

GO  TO  7b 

75  lja)C+l '  —5 1-ijM (.hB'3 ( PP)  +BOF'TrF)  .PP) 
V'.UC+l  >  =0. 

IUNC)  -OF'  U  (HC+ 1 ) 

vaio  =o. 

76  DO  77  1=1. MC 

77  Hi  I ) *8 < I +2 :  _  _ 

GO  Tn  4 

100  RETURN 
END 


(con t *) 

Co.vt'Aua  °'f  • 

foot  SaUfoUTme, 


« 


!(<•  32V9- 
10*  3**21  *9- 
l«.Tp  '«*,*5- 
10*  36 fe**5- 
10*3919*5- 
10*  Jf.Cr'S- 
I0*3I£2*S- 
t 0+361 o'S- 
iooz^^- 


Ul*332’*i1'- 

1 0+3286  •£- 

t0+3609•£- 

I0+36I6E- 

)0+3S»2*£- 

I0+3S92*E- 

I0+38E£*£- 

I0+366E*!- 

10+ 3216*1- 

10+3SS6*£- 

10+3256E- 

I0+33E6  1- 

I0+326£*£- 

I0»3<n*££- 

10+3622*1- 

16+326! *£- 

10+3601  *£- 

10+3600*1- 

1 0+3663  *2- 

10-3  222*2- 

!0+3£69*2- 

10+3866*2- 

!0+306£*2- 

10+3691  *2- 

10+3286*1- 

to*3ezz‘t- 

10+3655*1- 

10+3221*1- 

13+3260*1- 

00+3558*0- 

00+3185*2- 

00+3220*2- 

08+3562*2- 

00+3112*2- 

00+3261  *8- 

03+36E6*8- 

00+3309*8- 

00+  3669  *  0- 

00+3585*8- 

00+3926*8- 

00+3631*8- 

00+3968*2- 

08+3£26*2- 

00+3366*9- 

00+  3'J?2  *  9- 

00+3086*5- 

00+3525*6- 

00+3£6£*£- 

00+3958*1- 

20-3SIS*9 

00+ 3859 *Z 

00+3£l6*9 

i 0+330Z* I 

I0+369E* I 

00+3338*9 

00+3lr6*Z 

i 6-5505 *E 

00+3613*1- 

00+36!2*£- 

00+3181*6- 

00+3182*5- 

00+3396*5- 

08+3896*9- 

0O+3£08*9- 

00+3186*9- 

00+3210*2- 

00+3268*9- 

00+3519*9- 

00+3031 *9- 

00+3215*5- 

00+3989*6- 

90+ 3895* £- 

00+3001 *2- 

10-3068*1- 

80+36S£*2 

80+3125*5 

0O+3£22’2 

00+3305*5 

00+3910*2 

ie-38S£*8- 

O0+365PE- 

00+3110*6- 

00+3512*9- 

00+3006*2- 

00+3096*0- 

08+3651 *6- 
00+3^041*6- 
10+  3t£0*  1- 
10+3120*1- 
18+31011- 


20+3000*1 


10*3000*5 


*0 


10+3080*5 


20*3000*1- 


1«>04A 


£0+3080*6 
£0+JI96*C 
10+3I261 
£0+3280* £ 
£0+3268* £ 
E0+3£08*£ 
C0*3692*£ 
£8+ 3622*1 
£0+3589*1 

£0+3809*1 
£8+3295*1 
£8+322S*£ 
£0+3086 *£ 
£0+3866*£ 
£0+3606 *1 
£0+302£*£ 
£0+38££*E 
£8+3162*1 
£8+3252*£ 
£0+3212*1 
£0+3121  *£ 
£0*3££1 *1 
£8+3668*1 
£0+3550 *£ 
£0+3St0‘£ 
£8+3926*2 
£8+3916 *Z 
£8+3368*2 
£0+3850*2 
£0+3818*2 
£0*3622 *Z 
£0*36E2*2 
£0+3082*2 
£0+3199*2 
£0*31 29*2 
£0*3285*2 
£0*3265  2 
£0*3£0S*2 
£0*3696*2 
£8*3626*2 
£0*3581*2 
£0+3561*2 
£0+3981*2 
£0+3292*2 
£0*3322*2 
£0*3801*2 
£0*3861*2 
£8+3681*2 
£0*3028*2 
£8*3018*2 
£0*3166*1 
£0*3256*1 
£0+3216*1 
10+3E201 
£0+3££0*l 
£0+3662*1 
£0+3553*1 
£0+3512*1 
£0+3929*1 
£0+3919*1 
£0+3265*1 
£0+3855*1 
£0+3815*1 
£0+3626*1 
'■“+3616*  I 
10*9006*6 
£8+3191*  1 
£0+3IZ£"l 
£0+3282*1 
10*3262*1 
£0*3102*1 
£0*3631  *  I 
£0*3621*1 
£0*3580*1 
10*3560 * l 
£0+3900*1 
28+3299*6 
20+3122*6 
20+3628*8 
20+3586*0 
20+3160*8 
20+3269* 2 
Z0+3£0£*2 
20+3606*9 
20+3515*9 
20+3121*9 
20+3^22*5 
20+3£££*S 
20+3616*6 
20+3565*6 
20+3251*6 
20+ 1852*1 
20+3691 *1 
20*3826*2 
20*3925*2 
20*3281*2 
20+3002*  t 
20* 3661*1 
20*3406*1 


0302 


ViaC^ 


s  (i.-0/Ti 


(o')  D  DT 


Vmt*1 _ 


Cb^)  LOT 


Fiaufg  ^>G  ;  ~ Trans-(:o<<,'»'ti»i'5 

-ft,  dig.W  )|wplei»eK+^'«" 

o*f  Qnal ot\  Singly 

LC 


taa>  3'iniiH9H4 


10E+03  .  .87E+03  .  . ibE+04  .  .24E+04  .  .32E+04  .  .40E+04 

aqP4.f»7  17  F+Od  POF+Rd  . 2BE+04  .^hF+04 


IlHGUirUi'E  i  I'B  ■ 
14.61021 


iTUf.E 


Precision 


C.tk  L  %  &  9i  5  & 


IFkU 

G44.  7.0 

ltu>  Hl-I 

111!  4.1 

S'  bits 

G4G 

l.o 

|«*o  41.1 

xm  H.i 

"1  bit* 

G4G 

1.0 

IGOO  4H.8 

rwi  4.1 

G  b*ts 

GW 

G.l 

IS43  14  3 

xiS2  r\ 

S  bits 

GW 

“7.S 

mi  si.i 

2°}  4i  in 

4  b  rh 

Go 1 

3.9 

1177  4S.4 

2G35  U  2 

Ft  jure,  tT.  •  Res*Hs 

at  senSi't'i'v^  s4uJ>| 
-for  Str«j!y  4€»0Mn«M 
L.OX  tavtUUf 


Conclusion 


Three  basic  filter  structures  have  been  Investigated  for  implementing 
a  ten  pole  transfer  function  of  the  human  vocal  tract.  From  a  practical 
standpoint  only  the  lattice  filter  has  sensitivity  properties  good  enough 
to  warrant  serious  consideration. 

In  the  process  of  developing  the  synthesis  procedures  for  the 
specific  digital  filters,  several  more  general  computer  program  segments 
4ere  developed,  such  as  the  plotter,  the  complex  root  finder,  and  the 
?  algorithm  to  do  the  Cauer  expansion.  It  is  hoped  that  the  material  in 
this  paper  will  be  of  general  use  to  those  investigating  analog  sampled 
data  filters  in  the  future. 
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